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1 ABSTRACT 



2 The Cassini cameras have captured the opposition effect in Saturn's rings with a high 

3 radial resolution at phase angles down to 0.01° in the entire main ring system. We de- 

4 rive phase functions from 0.01° to 25° of phase angle in the Wide-Angle Camera clear 

5 filters with a central wavelength Aci=0.611 and phase functions from 0.001° to 25° of 

6 phase angle in the Narrow-Angle and Wide-Angle Cameras color filters (from the blue, 

7 Abi=0.451 fim to the near infrared Ait=0.752 /im). We characterize the morphology of the 

8 phase functions of different features in the main rings. We find that the shape of the phase 

9 function is accurately represented by a logarithmic model (Bobrov 1970, in Surfaces and 

10 Interiors of Planets and Satellites, Academic, edited by A. Dollfus). For practical pur- 

11 poses, we also parametrize the phase curves by a simple linear-by-part model (Lumme 

12 and Irvine 1976, Astronomical Journal, 81, p865), which provides three morphological 

13 parameters : the amplitude and the Half- Width at Half-Maximum (HWHM) of the surge, 

14 and the slope S of the linear-part of the phase function at larger phase angles. Our anal- 

15 ysis demonstrates that all of these parameters show trends with the optical depth of the 

16 rings. These trends imply that the optical depth is a key-element determining the physical 

17 properties which act on the opposition effect. Wavelength variations of the morphological 
IB parameters of the surge show important trends with the optical depth in the green filter 
19 (Agi.=0.568 /xm), which implies that grain size effects are maximum in this wavelength. 



20 Keywords: Saturn's rings; phase curves; opposition effect; coherent backscattering; shad- 

21 owing; shadow hiding 
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22 1 Introduction 



When one views tlie rings of Saturn with the Sun directly behind the observer, a phe- 
nomenon called the opposition effect can be seen. The opposition effect, also known as the 
opposition surge, is a sudden, nonlinear rise in brightness with decreasing phase (Sun- 
ring-observer) angle that occurs as the phase angle approach es zero. The opposition 



effect, which was observed for most Solar System bodies (e.g., 



Shkuratov et al.. 1999 



Verbiscer et al.. 2005 



Belskava and Shevchenko. 2000 



2007 



Helfenstein et al.. 1997 



Belskava et al.. 2003 



Simonelli and Buratti. 2004J: 



Rosenbush et al.. 20061 ) and (IRabinowitz et al.. 20071 ) was dis 



covered in t he course of Miiller' s long-term photomet ry of the Saturn system, beginning 



31 in 1878 (e.g. 



Miiller. 1885 



18931). 



Seehger (18841 



18871 ) inferred that the opposition effect 



32 was due to th e rings, since Jupiter did not show a comparable opposition brightening 



33 f lPoUack. 1975h . 



34 Most recent studies of the opposition effect in the Saturn's rings were based on ground 



35 based and spatial data which resolved t 



Poulet et al.. 2002 



le main rings (ILumme and Irvine. 1976 



Esposito et al.. 1979 



French et al.. 20071 ). Earth-based observations, though valuable, have 



37 a low to moderate spatial resolution (for the Hubble Space Telescope, 1 pixel at Sat- 



urn =285 km and the full-width at haf 
tion of Bubble's Camera =7 pixels, see 



maximum (FWH M) of the point-spread func- 



Cuzzi et al.. 20021 ). This indicates that the sig- 



nal is averaged over large regions in the rings. Unfortunately, the rings are highly het- 
erogeneous features th at may present rapid spatial variations of their optical properties 



( lEsposito et al.. 19871 ). So the interest of spacecraft observations is the ability to probe 
the signal in very narrow ring features (~40 km in the present paper), which may have 
much more uniform optical properties. This allows an easier study of the opposition effect 
which is an already complex phenomenon. 

Pioneer 11, Voyager 1 and Voyager 2 did not observe the rings at phase angles near 



zero. However, Cassini, which became an artificial satellite of Saturn in 2004, is the 
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first spacecraft to o 



49 flNelson et al. 2006 



jserve the opposition effect in the rings, with several instruments 



Altobelh et ah. 20071 ). including the Narrow Angle (NAC) and Wide 

50 Angle (WAC) cameras of the Imaging Science Subsystem (ISS). For our data set (Ta- 

51 blel), typically for the WAC 1 pixel=40 k m and the FWHM is 1.8 pixels, for the NAC 



52 1 pixel=5 km and the FWHM is 1.3 pixels (jPorco et al.. 20041 ). With Cassini images we 

53 are thus able to characterize how the surge varies throughout individual features in the 

54 ring system. 



The opposition effect is now known to be the combined effect of the coherent backscat- 
ter (at small phase angles) which is a constructive interference of photons in a medium 
of using grains of wavelength size and the shadow hiding (at larger phase angles) which 
consists of shadows cast by pa rticles on other particles that become invisible to the ob- 



59 server (IHelfenstein et al.. 19971). Some analytica. 



two effects ( IShkuratov et al.. 1999 



models have tried to combined these 



Hapke. 20021 ). but recent laboratory measurements 



and data analysis show that models' results could be compromised by some physical as- 
sumptions on the scatterer elements, the wave l ength-dependence and the the angular siz e 



63 of the source fIShepard and Helfenstein. 2007 



Shkuratov et al.. 2007 



Deau et al.. 20081 ). 



Two other effects, the external near-field effect and the single scattering internal-field 
coherence, were never been observed solely because they are intrinsically inseparable 
from the coherent backscatter effect and are for the moment only numerically simulated 



67 f lPetrova et al.. 2007 



Muinonen et al.. 20071 ) 



Both shadow-hiding and coherent backscattering are likely to play r oles in determining the 



69 shape of the phase curves of Solar System bodies at low phase angles (IHelfenstein et al.. 1997 : 



70 iHapke et al.. 19981 ). Shadow-hiding probably dominates at phase angles greater than a few 

71 degrees, while coherent backscattering takes effect at the very smallest phase angles. The 

72 shadow-hiding effect gives cl ues about the three-dimen sional structure of a layer of ring 



73 particles whic h, according to 



74 15m (see also 



Salo and Karjalainen. 2003 



Kawata and Irvine f 19741). could have a typical size of r 



French et al.. 20071 ). By contrast, the coherent 
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75 backscattering component sheds light on the nature of the ring particle surfaces at a scale 

76 not much larger than optical wavelengths. Indeed, previous photometric studies which 

77 investigate the opposition effect in the Saturn's rings determined ranges of typical size 



78 for coherent-backscattering effect of ab out d = 10 fim flMishchenko and Dlugach. 1992a 



Poulet et al.. 2002 



French et al.. 20071 ). Thus, the shape of the opposition effect provides 



80 cons traints on features v astly smaller than the camera's resolution. 

81 Like 



Poulet et al. (20021 ). we refer to particles as macroscopic and individual objects in 



82 the rings and grains as microscopic objects at the surface of rings' particles. Using this 

83 terminology, ring particles might be covered by regolithic grains. 



84 This is the first of a series of papers dealing with the opposition effect in Saturn's rings as 

85 seen by ISS/Cassini. In the present paper, we focus on the characterization of the morpho- 

86 logical modeling of the shape of opposition effect in the rings. Since the recent theories 

87 have difficulties to link the behavior of the opposition effect with the physic al proper- 



ties of the surface material ( Nelson et al.. 2000 



Shepard and Helfenstein. 20071 ). we think 



89 it is necessary to have an experimental approach, consisting in (1) deriving the shape's 

90 parameters and (2) trying to find correlations among themselves, and with the optical 

91 depth. In the second paper (Deau et al. in preparation), we will use recent analytical 

92 photometric models to derive some of the physical properties of the ring material. Indeed, 

93 to constrain completely the physical properties of the ring material, photometric and po- 

94 larimetric phase curves are needed, however, Cassi/ISS did not captured the opposition 

95 spot with its polarized filters. 

96 In section 2 of the present paper, we describe the ISS/Cassini data set and our procedure 

97 for extracting photometric data from images and fitting empirical models to the data. In 
93 sections, we characterize the morphology of the opposition surge at different locations in 

99 the main rings and focus our attention on cross-correlations among the morphological pa- 

100 rameters and correlations of these parameters with the optical depth and the wavelength. 

101 Finally, in section 4 we discuss these results and examine to which extents photometric 
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102 models may explain the morphological trends derived in section 2. 



103 2 Observations and reductions 



104 2. 1 The Cassini Imaging Data Set 



105 The Cassini ISS instrument is composed of two cameras, a wide angle camera (WAC) and 

106 a narrow angle camera (NAC) equipped with 1024x1024 CCD matrices. Both use a set 



of about twenty filters ranging from near-infrared to ultraviolet (jPorco et al.. 20041 ). 



Here, we focus only on filters in the optical domain, with first the blue, g reen, red and near 



109 infrared filters, that will be called hereafter COLOR filters. According to (iPorco et al.. 20041 ). 

no COLOR filters from NAC and WAC cameras do not have exactly the same central 

111 wavelength. For the blue filters, the NAC can have two different combinations of fil- 

112 ters (CL1,BL1) or (BL1,CL2) which lead to two central wavelengths : A^f^ = 0.440/im 

113 and A^f*^ = 0.451/im. In contrast, for the WAC, the blue (CL1,BL1) filter is charac- 

114 terized by a wavelength of about A^"^*^ = 0.460/xm. Because the spectral width of these 

115 filters is about ±0.050/im, we consider all images from blue filters of the NAC and the 

116 WAC as a consistent whole. The green filter (CL1,GRN) has almost the same spectral 

117 characteristics for the two cameras since the central wavelength of Ag^jf = 0.568/im for 
iiB the NAC and A^'^,'' = 0.567/im for the WAC. The red filter for the NAC corresponds to 

119 the filter (RED,CL2) with a central wavelength of A^^^^"" = 0.650/im. For the WAC, the 

120 combination (CL1,RED) is at a central wavelength of X^^^ = 0.649yum, which does not 

121 change to the that of the NAC at one nanometer. Finally, the filter (CL1,IR1) in the 

122 near infrared shows a difference of 10 nanometers between the central wavelength of the 

123 NAC (Af,*'' = 0.752/im) and that of the WAC (A**^ = 0.742/im). In summary, with this 

124 moderately high spectral resolution, the combination of images that are not coming ex- 

125 actly from the same filters are responsible for a shift of the central wavelength from 0.001 
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126 to 0.010 micron. Because each filter has a non neghgible spectral range (±0.050/xm), we 

127 consider without distinction the COLOR filters from the NAG and WAG cameras and in 
12B the rest of the paper, we take as a reference the following wavelengths (Abi = 0.45 1/im; 

129 Agrn = 0.568/xm; Ared = 0.650/im and Air = 0.752yum) to designate each GOLOR phase 

130 curve. 

131 The other part of images discussed in the present paper has been taken in GLEAR mode 

132 of the WAG, designating the absence of filters resulting in a spectral bandwidth spanning 

133 from 0.20 to 1.10 /im (the central wavelength is Ad = 0.611/im and the spectral range is 

134 ±0.450/im). 

135 For the present paper, the selected observation campaign of the opposition effect on the 

136 Saturn's rings is divided in four sequences from May 2005 to July 2006 (Table 1). 

137 [Table 1] 

13B These observations were conduced conjointly by the VIMS, GIRS and ISS instruments. 

139 Three of ISS sequences were obtained using the WAG only and the last one using the 

140 BOTSIM mode (using both cameras simultaneously). Depending on the sequence, different 

141 filters were used and different spatial and angular resolutions were achieved. In the present 

142 paper, the first of a series aimed at a detailed study of the opposition effect, we will focus 

143 on sequences obtained using the WAG camera in GLEAR filter mode and sequences using 

144 NAG and WAG cameras in GOLOR filter mode. Table 1. Images of the B and G ring are 

145 shown in Fig. 1 and 2 to illustrate the quality of the data set. 

146 [Fig. 1, Fig. 2] 

147 In the GLEAR filter mode, the wide angle camera captured rings at zero phase angle two 
14B times in 2005, on June 7 and June 26 (they will be designated as June 7 and June 26 

149 sequences in the rest of the paper, for practical purpose). 

150 In the June 7 sequence, the WAG images have an average resolution of ~ 44 km per 

151 pixel; the rings are observed in reflexion. In each individual image, the phase angle varies 
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152 by about 3 degrees. In the majority of the set, the opposition point at zero phase an- 

153 gle is visible with a very good samphng : about 0.005 degree per pixeL This set covers 

154 the full ring system by tracking the opposition spot. Most of the set has an excellent 

155 photometric quality, however in dark regions the strong Saturnshine produces ghost im- 

156 ages of the secondary mirror perturbing the signal significantly. The amplitude of this 

157 artifact is quite c onstant, about 0.05 in I/F (I/F is the phase-corrected reflectivity, see 



158 iPorco et al.. 20041 ). which does not significantly affect the photometry of the bright re- 

159 gions as A and B rings. However it is very troublesome for dark regions like the D ring, the 

160 C ring and the Cassini Division. This is why images from the June 7 sequence have not 

161 been considered f or these two regions. The same instrument artifact has been observed 



162 and discussed in (iHedman et al.. 2007a[ ). 

163 The June 26 sequence has similar characteristics, see Table 1, with a better radial reso- 

164 lution of about 30 km per pixel and spans the B and C rings. The presence of the ghost 

165 signal in the C ring is not clear. A detailed photometric analysis in these images (cf. 

166 section 3.1.1) shows that the derived phase function is consistent with an unperturbed 

167 signal, very differently from images in which the ghost were clearly identified. It is why 
16B the images of the C ring taken on June 26 have been considered for the C ring phase 

169 function. 

170 The July 23 sequence has the best radial resolution (~13 km.pixel"^) and spans all A, B, 

171 C rings and the Cassini Division (the radial location of the opposition spot in each image 

172 is represented in the Fig. 3). In this set, the ghost artifact is absent owing to a much larger 

173 angular separation with the bright Saturn globe. Some images of this sequence provide 

174 also the larger phase angles. 

175 [Fig. 3] 

176 In COLOR filters mode, the ISS instrument uses the BOTSIM mode in the May 20 

177 sequence, using both cameras simultaneously for a ride in A and B rings, Cassini Division 
17B and the outer C ring (Fig. 3). However, because the boresight of the NAC and the NAC 
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179 are not exactly the same, NAC images did not represent a zoom of the opposition spot in 

180 WAC images. Thus, WAC and NAC images captured at the same time have opposition 

181 spots that are not exactly at the same location in the rings. The radial resolution in NAC 

182 images is 4.6 km.pixel"^ and spans from 44 to 66 km.pixel"^ for WAC images in COLOR 

183 filters. The other sequences {Dec 31, Feb 20 and Apr 25), did not have the opposition 

184 spot but provided the larger phase angles. 



185 2.2 Data reduction 



186 2.2.1 Calibration 



187 Raw images are calibrated first with a standard pipeline described in (jPorco et al.. 20041 ) 

188 and called CISSCAL, using the 3.4 version. The DN values are converted into the I/F 

189 units. This ratio is dimensionless, with / is the specific intensity measured by Cassini, and 

190 irF is the incident flux from the Sun. So it is a measure of the local reflectivity under the 

191 current observing geometry. 

192 Cassini images are not directly exploitable to study the photometric behavior of rings. 

193 In order to do this, the relevant data we need are the so-called phase function which 

194 is linked to the I/F ratio as a function of the phase angle, and corrected from effects 

195 of observation geometry. A full procedure has been designed to reconstruct the phase 

196 function from different images with different observation geometry and resolutions and is 

197 detailed below. 



198 2.2.2 Extraction procedure 



199 The first step is to reproject the images in a (Radius, Longitude) frame, in which features 

200 at a same radius from Saturn are horizontally aligned. This procedure critically depends 

201 on the quality on the nav igation. When possib le, the edges of the A, B, C rings as well as 



202 ring features reported in (lEsposito et al.. 19871 ) were used as fiducial references. Distances 
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203 reported in (lEsposito et al.. 19871 ) were corrected according to (iNicholson et al.. 19901 ). 

204 The resulting navigation error on WAC images varies from 1 to 1.5 pixels from the center 

205 to the edge of the image. 

206 The 1/F in a ring at constant distance from Saturn is obtained by extracting the data 

207 on a line of constant radius {i.e. an horizontal line) in the reprojected map of calibrated 
20B brightness. Other geometrical parameters are extracted in the same way : phase angle a, 
209 cosine of incidence angle /iq, cosine of the emission angle fi, optical depth of the r ings ob- 



210 tained from the PPS Voyager instrument (A = 0.260 p^m). (lEsposito et al.. 19871 ): radius 



Nicholson et al. f 19901 ). 



211 scale has been corrected with the procedure of 

212 This procedure works very well for structures with width larger than the navigation ac- 

213 curacy. However radial structures are visible at all scales down to one pixel, (Fig. 2). For 

214 structures radially smaller than 5 pixels, navigation errors make false the extraction along 

215 a hne of constant radius close to the edges of the image, inducing accidental extraction 

216 in nearby different features. To overcome this problem and to ensure that we always ex- 

217 tract the same ring feature, we developed a ring tracking technique using a basic pattern 

218 recognition algorithm to follow a single feature. Extensive visual check of the result shows 

219 the method is reliable down to 1 pixel of radial width. 



220 2.2.3 Construction of the phase function 



221 The ultimate information we need is the phase function of individual ring particles to 

222 characterize their surface properties. Unfortunately, the sign al from an individu al particle 



223 is heavily altered because of the finite thickness of the rings ( iPorco et al.. 20071 ). Also, the 

224 signal's intensity depends on observation angle with respect to the ring's normal. Inverting 

225 such complex collective photometric effects would require the use a detailed light scattering 

226 code with many assumpt ions concerning th e photometric prope rt ies of particles. Such code 



227 has 

228 and 



jeen developed by 



Salo ri992 



19951) 



Richardson (19941 ) 



Porco et al. fl999 



20071) 



French et al. (20071 ). However, for our present purposes, they cannot used to de- 



ll 



229 rive the phase functions in hundreds of different regions as we wish to do here. Conse- 

230 quently, as a first app roximation, we use in the present paper the classical approach of 



Chandrasekhar (1960l ) linking the I/F to the phase function with the following assump- 



232 tions : homogeneous layer of particles and singl e scattering. This 



after assumption may be 



233 justified for phase angles smaller than ~30° (ICuzzi et al.. 20021 ). In reflexion, the phase 

234 function vuq ■ -P(tt) is derived from the solution to the radiative transfer equation (the 

235 designation phase function is not strictly accurate since what we really determine is the 



236 product of the single scattering albedo, vjq, times the particle phase function, -P(a)) : 



zon ■ P(a) 



I 4(/x + /ip) 

— X X 

F /io 



M MO 



with r, /i, fio, a standing for : the normal optical depth, cosine of emission angle, cosine of 
incidence angle and phase angle respectively. In order to allow future comparisons with 
detailed numerical models, numerous analytical parametrization of the observations are 
provided by morphological model in section 2.3. As one can see, the value of the optical 
depth is necessary to derive the phase function. Preceding works has shown that the 
exponential factor can be neglected in first approximation for Earth-based observations 



243 fICuzzi et al. 2002 



Poulet et al.. 20021 ). However, we noticed that we obtained much more 



244 coherent results when taking into account the exponential factor, when comparing results 



Chandrasekhar fl960 


), as previous photometric studies based on 


of Saturn's rings 


(Dovle et al.. 1989: 


Cooke. 1991: 


Ferrari. 1992: 


Dones et al.. 1993 


. 



spacecraft observations 



Showalter et al.. 1992 



249 -P(«) obtained with Eq. 1 is the particle's disk integrated phase function which determines 

250 the angular distribution of single scattered radiation from the body as a whole. The 

251 phase function is normalized over the solid angle Q to the single scattering albedo : 

252 wq = ^ J P{a) dfl. To derive the albedo, the full phase curve, from to 180 degrees of 
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253 phase angle, must be known. In this study, we have restricted our data (0<a;<25 degrees), 

254 so that we may avoid separating the albedo from the phase function. The derivation of 

255 the albedo will be presented in the second paper which relates ISS data from to 180 

256 degrees. 



257 2.2.4 Finite size of the Sun 

258 Generally speaking, all phase curves present a bright opposition surge below 1° and a 

259 slope decreasing linearly for a >1°. Whereas WAG images go down to ~0.01° of phase 

260 angle (Fig. 1), NAG images taken in GOLOR filters capture the opposition spot at better 

261 angular resolution (Fig. 2), the resulting phase functions from NAG images go down to 

262 ~0.001° (Fig. 4). This is the first time the opposition spot is imaged at such fine scale. 

263 Note that we define the phase angle as the angle between the vector pointing to the Sun's 

264 center, and to the spacecraft, from the observed point. So from a strictly mathematical 

265 point of view there is no lower limit to the phase angle value despite the source's finite 

266 size. For an extended illumination source, the phase curve should be the integral of a point 

267 source phase function over the Sun angular radius. Thus on a more physical point of view 
26B i t is not possible to observe the phase func tion at angle below the Sun's angular radius 



269 ( Kawata and Irvine. 1974 



Shkuratov. 19911 ). 



270 Fig. 4 demonstrates that the opposition surge flattens (in all rings and at all wavelengths) 

271 at phase angles below 0.029°, in good agreement with the Sun's angular radius at the 

272 date of observations (0.0291 degree, given by a^^i^ = arcsin „ — where Tq = 6.96 x 

-f^Q — Saturn 

273 1 0^ km is the radius of th e Sun and -Ro-satum is the heliocentric distance of Saturn 



Murray and Dermott. 200 



[Fig. 4] 
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276 2. 3 Data fit with Morphological models 



277 The purpose of the present paper is to provide an accurate description of the morpho- 

27B logical behavior of the observed phase curves. This is the very first step prior to any 

279 attempt of further modeling, which will be the subject of a forthcoming paper. As a con- 

280 sequence, special care has been given to parametrize the phase functions conveniently. In 

281 addition morphological parametrization is necessary to compare efficiently hundreds of 

282 phase curves at different locations in the rings and derive statistical trend as will be done 

283 in sections. 

284 Several morphological models have been us ed in the past to describe the shape of the 



285 phase functions : the log arithmic model of 



Bobrov (1970l ). t he linear-by-part model o f 



Kaasalainen et al. (200ll ). 



286 iLumme and Irvine (19761 ) and the linear-exponential model of 

287 The specific properties of these three models make them adapted for different and com- 

288 plementary purposes. The logarithmic model is interesting for direct comparisons with 

289 numerical models, the linear-by-part model is convenient to describe the shape in an in- 

290 tuitive way, and finally the linear-exponential model is adapted for comparison with other 

291 studies previously published. 



292 2.3.1 The logarithmic model 



As 



Bobrov mm 



Lumme and Irvine (19761 ) and 



Esposito et al. (19791 ). we remark that 



294 a logarithmic model describes very well the solar phase curves of the Saturn's rings. It 

295 depends on two parameters (oq and ai). This first morphological model has the following 

296 form : 



WqP^o) = aQ + ai X log(a) (2) 
297 In general, this model is the best morphological fit to the data. It is reasonably accurate 
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29B down to 0.029° of phase angle. This is due to the finite angular size of the Sun which 
299 fiattens data below 0.029°, whereas the logarithmic function continues increasing (see 



300 section 2.2.4 and Fig. 4a' 

301 this was also noticed by 



. For large phase ang 



es, the fit is satisfactory up to a ^15 degrees, 



Altobelli et al. (20071 ) who fitted temperature phase curves of the 



302 C ring with CIRS/Cassini. 

303 Whereas the meaning of and ai are not easily interpretable in term of shape, to allow 

304 an easy comparison with future numerical simulations the values of these two parameters 

305 are reported in Table 1 of Electronic supplement material. 

306 2.3.2 The linear- exponential model 

307 This model describes the shape of the phase function as a combination of an expo- 
30B nential peak and a linear part. Its main interest is that it has been used in previ- 



309 ous work for the study o : 



310 (IKaasalainen et al.. 2001 



the backscattering of Solar System's icy satellites and rings 



Poulet et al.. 2002 



French et al.. 20071 ). We give the details of 



311 this model : the 4 parameters are the intensity of the peak Jp, the intensity of the back- 

312 ground /{,, the absolute slope of the linear part Ig and the width of the exponential w, 

313 such that the phase function is represented by : 



wqP^o) = h + Is - a + Ip - e 



(3) 



314 With these 4 parameters, we characterize the shape of the phase function by introducing 

315 three morphological parameters : A, HWHM and S designating the amplitude of the surge, 

316 the half-width at half-maximum of the surge and the absolute slope at large phase angles 

317 respectively, so that : 



A 



Ip + h 



HWHM = 2 • In 2m; and S = -J, 



(4) 
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As iFrench et al. (20071 ). we noticed that this model did not fit well the phase curves, in 
particular, the derived model parameters appeared to depend substantially on the phase 
angle coverage (see section 2.3.4), preventing a robust com parison of data. This is d ue to 



321 the fact that we did not use the converging procedure of 



Kaasalainen et al. (200ll ) but 



322 rather the common downhill minimization technique, as done by t he previous users of 



323 the linear-exponential model (jPoulet et al.. 2002 ; 



French et al.. 20071 ). Moreover, it seems 



that the angular scale at which the phase function is observed may strongly influence 
the model parameters. So to avoid this problem, we used a much simpler morphological 
model that appeared muc h more robust for the com parison of heterogeneous data set : 



327 the linear-by-part model of 



32B of the data with the linear-by-part model of 



Lumme and Irvine f 19761). also we found a much better match 



Lumme and Irvine (19761 ) (see below) 



329 However, the linear-exponential was consid ered in the present study only to unde rstand 



330 variations of A, HWHM and S between the (jPoulet et al.. 20021 : 



French et al.. 20071 ) phase 



331 curves and our phase curves (see section 2.3.4 and section 3.2.3). 



332 2.3.3 The linear-by-part model 

333 For an intuitive description of the main features of the phase curves, the linear-by-part 

334 model is the most convenient one. It is constituted of two linear functions fitting both the 

335 surge at small phase angle (a < ai) and the linear regime at higher phase angle (a > 0^2) : 

zuoP{a < ai) = —Aq x a + BozuoP^a > 02) = —Ai x a + Bi (5) 



Lumme and Irvine (19761 ) and 



Esposito et al. (19791 ) use ai=0.27° and ^2=1. 5°. However, 



337 we encountered difficulties with the value of ai. By testing several values of ai, it appears 

33B that for our data set, values of ai less than 0.3° yield a general overestimation of Aq, 

339 especially in the C ring and values of ai greater than 0.3° yield an underestimation of 

340 Aq only in the B ring. Consequently we found the our data were better reproduced using 

341 ai = 0.3° which is now adopted in the rest of the paper. 
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342 With this four outputs : the two absolute slopes Aq and Ai and the two y-intercepts Bq and 

343 Bi, the shape of the curve is characterized by the morphological parameters A, HWHM 

344 and S are then defined by : 



A=:^ HWHM= and S = A, (6) 

345 The purpose of this model is not, of course, an accurate description of the data but rather 

346 a convenient description of the main trends of the phase curve. 



347 2.3.4 Stability of the morphological parameters 



Poulet et al. f200i ) and 



French et al. f2007h 



348 In order to compare properly our results with those ofO 

349 for which the Saturn's rings phase curves did not have the same phase angle coverage, we 

350 have tested the influen ce of the portion 0.0 5°-0.4° and 6°-25° on the converging solution. 



351 In deed, phase curves o f 



of 



Poulet et al. (20021 ) and 



Poulet et al. f2002[) do not have data under 0.3°, and both studies 



French et al. (20071 ) do not have data over 6.5° (the maximum 



Ft 

arcsin 0-^=^^^ < g.50 where i?0_Earth is 

-'^© — Saturn 



353 phase angle reached from Earth is about aQ^nax 

354 the heliocentric distance of the Earth and Rg^-fj^tum is the heliocentric distance of Saturn 



355 computed with orbits of 



Murray and Dermott. 2000l ). 



356 With two typical Saturn's rings phase curves of ISS (Fig. 5), we have removed data by sec- 

357 tion of 0.1° and fit the pseudo-incomplete phase curve with the linear-exponential model, 
35B which provides the new solution, designated by ^remove, HWHMremove and Sj-cmovc- The 

359 initial solution found for fuller phase function (0.01°-25°) is called ^optimal, HWHMoptimai 

360 and S'optimai- Both solutions are ob tained with a down hill minimization techni que, to re- 



361 produce the fi tting method used by 



362 to recall that 



Poulet et al. f2002h and 



Kaasalainen et al. f200ll') 



363 however because (e.g. 



French et al. (20071 1 . It is wise 



proposed a convergi ng procedure more accurate. 



Poulet et al.. 2002 



French et al.. 20071 ) did not used it, and because 



364 we want to reproduce not the best converging solutions but the deviations lead by the fits 
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365 driven by iPoulet et al. (20021 ) and lFrench et al. (20071 ), that is why we use the downhill 

366 method and not the probability distribution method of iKaasalainen et al. (200ll ). 



367 [Fig. 5, Fig. 6] 

36B In Fig. 6, we represent the ratio of the morphological parameters of the incomplete phase 

369 curve over the morphological parameters of the fuller phase curve, called v4rcmovcMoptimab 

370 HWHMrcmovc/HWHMoptimai and S'rcmovc/'S'optimai- We observe a slight underestimation of 

371 A, a strong overestimation of HWHM and a moderate underestimation of S. 

372 The deviation of the optimal value is quite weak for A (ArcmovoMoptimai ~0.96 at a cut- 

373 off of 0.3°) but its variation depends on the morphology of the surge. For the typical 

374 B ring phase curve (which has a narrower peak), we note a slight decrease of the ratio 

375 ArcmovcMoptimai and for the C ring, we note a slight increase of the ratio. However the 

376 both ratios lead to an underestimation, which means that incomplete data fitted by the 

377 morphological model will have a smaller amplitude that data which cover the full surge. 
37B HWHM shows the strongest deviation of Aromove- Indeed, HWHM is overestimated for 

379 the typical B ring phase curve (HWHMrcmovc/HWHMoptimai ~1.4 at a cutoff of 0.3°) but 

380 is strongly overestimated for the typical C ring (HWHMrcmovc/HWHMoptimai ~2.0 at a 

381 cutoff of 0.3°) which has a wider peak. 

382 For the slope, we notice an overestimation in the order of S'removc/'S'optimai ~1.3 at a cutoff 

383 of 7° for the B ring phase curve and of ~2.5 for the C ring at the same cutoff. This means 

384 that the slope of the linear part is stabilized at roughly 15°. 



385 These comparisons are important when we compare morphological trends found by 



386 in section 3.1.2 and 



French et al. (20071 ). in section 3.2.3. 



Poulet et al. f2002h 



387 2.3.5 Linking morphological parameters with the physical parameters of the models 

388 The use of a simple morphological model is generally not adapted to derive the physical 

389 properties of the medium. However, the theories developed for the coherent backscatter- 
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390 ing and the s. 

391 phase curve 



ladow hiding effects deduce their 



Hapke. 1986 



Mishchenko and Dlugach. 1992a 



pr operties by parameterizing the opposition 



Mishchenko. 1992 



Shkuratov et ah. 1999 



20021 ). Recent laboratory experiments raise some doubts on the meaning of 



393 the physical parameter s of these models and their corr e lation with the real phy sical prop- 



394 erties of the medium (IShepard and Helfenstein. 2007 



Shkuratov et al.. 20071 ). However 



395 because there is no better modelization, we connect the morphological parameters A, 

396 HWHM and S with the physical characteristics of the medium derived from these models. 



A is the amplitude of the opposition peak and describes the behavior of the phase func- 
tion at the smallest phase angles (a<2°). What we know about the opposition ef- 
fect is that it occurs at the smallest phase angles and acts on the multiple scattered 
light in the regolith on the surface of the particles : the underlying phenomenon is the 
coherent backscattering effect. The coherent-backscattering effect increases in bright- 
ness by almost a f actor two, while using grains size smaller than the wave length of 



the incident light, flMishchenko and Dlugach. 1992b 



Shkuratov et al.. 19991 ). In con- 



trast, the second phenomenon of the opposition effect, the shadow hiding effect, is 
known to produce a wide peak from to 2 degrees, a nd to decrease the brightness 



406 up to 20 degrees (IHapke. 1986 



Stankevich et al.. 19991 ). The combination of the two 



effects at very low phase angle is still a matter of debate and tod ay two theories dis 



408 agree in order to explain the pea. 



i of the opposition. The theory of 



Mishchenko fl992h : 



Mishchenko and Dlugach (1992bl ) assumes t hat the opposi tion peak is a pure coherent 



Hapke (20021 ) shows that the opposition 



410 backscattering effect whereas the theory of 

411 peak results from a coupling of coherent-backscattering and shadow hiding, even at low 

412 phase angles. This coupling should be due to the fact that the coherent backscatter- 

413 ing could act on both multiple and single sc attered light whereas the shadow hiding is 



a single scattered light effect (IHapke. 20021 ). Thus, this theoretical model defines two 



amplitudes as output parameters : the coherent backscatter amplitude Bqo and the 
shadow hiding amplitude Bso- As a consequence, using this theory, it does not seem to 
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be possible to ascribe the morphological parameter A solely to the coherent backscat- 
tering. 



419 For most laboratory measurements (jShkuratov et al.. 1999 



Nelson et al.. 200(1 ) . the 



amplitude of the opposition peak is a function of grain size in such way that A decreases 
with increasing grain size. This anti-correlation finds a natural explanation by the fact 



that macroscopic irregularities (>A) create 
ularities (< A). Laboratory measurements of 



ess coherent effects than microscopic irreg- 



opp osition surge increa ses when irregularities are small. 



Kaasalainen f2003l) also confirmed that the 



and 



Mishchenko and Dlugach (1992bl ) 



Mishchenko (19921 ) underline the fact A is linked to the i ntensity of the backg round 



lb which is a decreasing function of increasing absorption (ILumme et al.. 1990l ). thus 



A must increase with increasing absorpti on or decreasing a lbedo Wq, which was also 



Kaasalainen (20031 ) who found that the peak 



confirmed by laboratory measurements of 
decreases with increasing sample albedo. 
HWHM , the half-width at half-maximum, is generally associated to the coherent backscat- 
ter effe ct. It can be related to the grain size, index of refraction and packin g density of re- 



432 gohth flMishchenko. 1992 



Mishchenko and Dlugach. 1992a 



Hapke. 2002h . The HWHM 



is maximum for a effective grain size near A/2 and increases when the regolith grains fill- 
ing factor / increases. For high values of /, the H WHM shifts to wards the greater grain 
size. However, as for the amplitude, the model of lHapke (20021 ) defines two HWHMs : 



th e coherent backscat ter angular width {he)-, which is defined similarly that in the model 



of iMishchenko (19921 ). and the shadow hiding angular width {h^). This reinforces the 
idea that the observed surge results from a coupling between the coherent backscattering 
and the shadow hiding. 
S , the slope of the linear part, seems to be due only to the shadow hiding effect: the 
interferences caused b y coherent backscattering effect seem not to be very significant 



442 at larger phase angles (IMishchenko et al.. 2006 



Hapke. 20021 ). This means the shadow 



hiding light is not affected by the coherent backscattering at larger phase angles. 
Also, according to recent analytical and numerical models, the shadow hiding acts solely 
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on the linear part of the phase function (IStankevich et al.. 1999 



Shkuratov et al.. 19991 ). 



We define here the absolute slope, whereas it is naturally a negative parameter, be- 



cause the phase function a . 



( iKawata and Irvine. 1974 



ways decreases in brightness from 10° to 40 ° of phase angles 



Stankevich et al.. 1999 



French et al.. 20071 ) 



he slope depends on the volume filling factor D and the optical depth of the slab 



( IKawata and Irvine. 1974 



Stankevich et al.. 19991 ). Two models has been developed for 



different regimes of the particle volume density D and predict opposite behavior of S 
as the function of D and r. 
The shadowing model of 



particle shadow hiding by 



rvine (19661) and 



Kawata and Irvine (1974J ). also called inter- 



scopic par ticles (r=15 m, see 



Goldreich and Tremaine (19781) . consists in shadows of macro- 



Kawata and Irvine. 19741 ) in a particulate medium, such 



as 8D<Cl (jlrvine. 1966( ). The sm aller the volume density D , the steeper the phase func- 



model exists (lEsposito. 1979 



tion for incre asing phase angle (IKawata and Irvine. 19741). Other refinements of this 



Hapke. 1986 



Stankevich et al.. 1999 



French et al.. 2007h 



but lead to the same results : the opposition peak due to the shadowing sharpens and 
the absolute slope increases with decreasing packing density (or filling factor). 



Another model exists, the intra-pa rticle shadow hiding, which is y a. 



tide volume density accor ding to 



Goldreich and Tremaine (1978[ ): 



id for higher par - 



Muinonen (19941 ). 



Buratti and Veverka (19851 ) underlined the fact that the mutual shadowing among re- 
golithic grains could be suited for understanding the textural properties of the regolith. 
As a consequence, this mechanism operates at the surface of ring material : e.g. in the 
regolith layer. 

However, other physical parameters need to be ta ken into account. In the analytical 



inter-particle shadow hiding model of 



Hapke (20021 ). the slope can be linked to the an- 



gular width of the shadow hiding hsh- Thus a normalized and absolute slope would 



be S 



1 



Q t(A f)f in(i-D) ' "with D the volume filling factor, Qcxt(A,'r) the ex- 
tinction coefficient and f the mean radius of particles. Then, the slope should depend 
on the wavelength and on the particle size. Interestingly, laboratory experiments of 
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Kaasalainen (20031 1 showed that the slope increases when the sample's size increases. 
However, the microscopic and/or macroscopic roughness of the me dium need also to b e 



475 taken into account in t h e shadow hiding mo dels, as underlined by 



Shkuratov et al. f 19991 ): 



Kaasalainen f2003[) and 



477 example, laboratory experiments of 



Hapke (1984 



19861); 



Shepard and Helfenstein (20071 ). For 



Kaasalainen (20031 ) showed that the slope increases 



when the sample's roughness increases. 

Then laboratory measurements and recent theoretical models could significantly in- 
crease the number of physical parameters on which S depends. 



481 From the above arguments, the HWHM and the amplitude A are governed by both co- 

482 herent backscatter and shadow hiding effects whereas the slope S provides information of 

483 the shadow hiding effect solely. 



484 3 Shape of the phase curves at opposition 



485 3. 1 The opposition effect in CLEAR filters 



486 Due to the automation of extraction and fitting procedures (cf section 2.2) and due to 

487 the high images resolution, phase functions were extracted in as many as 211 different 

488 locations, in the D, C, B rings, Cassini Division and A ring (in increasing distance from 

489 Saturn, left column of Table 1 of Electronic supplementary material). In this section, we 

490 first present the typical behavior of some selected phase curves in different regions of the 

491 main rings (section3.1.1), then we discuss similarities and differences and what are the 

492 general trends from sections 3. 1.2 to 3.1.4. 
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493 3.1.1 Overview of the CLEAR phase curves 



494 Examples of phase curves in various ring regions are presented in Fig. 7. Each pair of 

495 graphs show on the left side a zoom from 0.01 to 2.5 degrees and on the right side, 

496 the fuller phase curve from to 25 degrees. These curves were obtained by combining 

497 several WAC images with a large distribution of viewing geometries (Table 1) each curve 
49B is built from the merging of 10 to 70 different images with various values of emission angle, 
499 incidence angle, phase angle etc. The dispersion of points is not due to the measurement 



Chandrasekhar (19601 ) inversion formula 



500 uncertainty but reflects mainly the limits of the 

501 that was used to extract the phase function from the measured values of I/F. It seems 

502 that some important physics may be missing (like multiple scattering), that could explain 

503 the scattering of points. Curves for the A and D rings (Fig. 7) are incomplete between 3 

504 to 20° which is due to removal of images because of an artifact (cf section 2.1). 



[Fig. 7] 



506 Whereas the general shape is similar from one ring to another (Fig. 7), some details in the 

507 shape may vary significantly. First pair of graphs in Fig. 7 shows the phase curve derived 
50B in the D ring from images of June 26. Due to short exposure time (10 ms), the D ring 

509 ringlets are too faint to be detected. In images of June 26, a bright spot is visible from 

510 67 000 km to the inner boundary of the C ring : this corresponds to the expected location o f 

511 the background sheet of material constituting part of the D ring (IHedman et al.. 2007al ). 

512 However camera artifact may be visible in such dim regions of the image but could not 

513 be clearly identified here. So the fact that a strong increase of brightness at the expected 

514 location of the opposition and the coherent variation of the signal with observing geometry 

515 between different images suggest that we indeed see the opposition effect in the D ring. 

516 However some doubts still remain. From 0.5 to 2 degrees, the curve is similar to other 

517 rings. Below 0.5 degree an exponential surge and a flattening at zero degree distinguish 

518 this phase function from the other ones. Does it reflect optical properties of D ring dust ? 
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519 
520 
521 



Is it an artifact ? This plateau below 0.5 degree is much too large to be explained by 
the finite angular radius of the Sun (0.025 degree). Because of these uncertainties, at this 
point it is speculative to interpret this specific behavior as real. 



522 For the C ring (Fig. 7), the shape of the phase function is well sampled below 2 degrees. 

523 It is comparatively wider than in dense A and B rings. More precisely, HWHM of the 

524 opposition surge is wider for the C ring (HWHM=0.26°) than for the A and B ring 

525 (HWHM~0.20°). This could be also interpreted as a steeper slope of the linear regime of 

526 the phase function for a > 2°. Wavy features between 5 to 25 degrees can be attributed 

527 to images artifact. Their amplitude is about 15% of the total signal of the C ring. 

52B The B ring opposition surge has the smallest a mplitude of all rings ( A=1.25 in Fig. 7). This 



529 was already underlined by the rec ent s tudy of 



530 servations of 



Esposito et al. (19791 ) and 



French et al. (20071) whereas previous ob- 



Poulet et al. (20021 ) give the exact opposite trend. 



531 Ho wever, their result c ould have a bi as due to the lack of d ata belo w 0.3° of phase angle 



for 



Poulet et al. (20021 ) and 0.1° for 



Esposito et al. (19791 ) whereas 



French et al. (20071 ) 



533 had values much smaller (the smal 

534 range from 0.0037 to 0.0132 degrees 



est phase angles at which HST observed the B ring 



French et al.. 20071 ). The B ring has also the steepest 



535 slope (S=0.105 WqP. deg ^) in the linear regime explaining why the opposition spot is so 

536 contrasty in the ISS images. 



537 The Cassini Division has a similar amplitude and width (A~1.47 and HWHM~0.28° in 
53B Fig. 7) to the C ring (A~1.45 and HWHM~0.26° in Fig. 7). Their slope of the linear regime 
539 are also similar (S=0.033 wqP -deg"^ for the Cassi ni Division and S=0 .030 ujqP .deg~^ for 



Poulet et al. (20021 ) and are suggestive 



540 the C ring). These similarities were first noticed by 

541 of a strong dependence of the opposition effect on the optical depth (we will come back 

542 to this in section 3.1.4). 



543 An example of the A ring phase function is given in the last pair of graphs in Fig. 7. 

544 At first sight, the signal appears much more disturbed than in other rings : specifically. 
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545 pieces of phase functions extracted from different images show a wide dispersion in this 

546 graph, whereas signal from an individual image has a very low dispersion. The origin of 

547 this is not clear and may be due to the artifact reported in section 2. The dispersion 
54B of the data is about 15% of the signal whereas the camera artifact should represent at 

549 most 5% of the signal only (estimated on the background). It may be possible that the 

550 dis persion may be also du e to an intrinsic photometric effect which is not corrected by 



551 the IChandrasekhar (1960l ) single scattering model (equation(l)). Indeed, the A ring has 

552 an intermediate optical depth ~0.5 so that it neither appears as a solid surface (like the 

553 B ring) or as a dilute system (like the C ring). Here, we are in an intermediate regime 

554 where many collective effects may influence strongly the apparent phase function (multiple 

555 scattering, gravitational w akes, density waves , etc.). A sophisticated model is required here 



556 to investigate such effect (jPorco et al.. 20071 ). However the general trends are quite clear 



557 and the A ring phase curve has a larger peak amplitude (A=1.39 in Fig. 7) than the B ring's 

55B (A=1.25 in Fig. 7). We see in addition that the slope of the linear regime is shallower in 

559 the A ring (S=0.078 zuoP.deg~^ in Fig. 7) than in the B ring (S=0.105 zuoP.deg~^ in 

560 Fig. 7) but steeper than in less dense rings (S~0.03 rooP.deg"^ in Fig. 7). Thus the phase 

561 curve at opposition in the A ring is somewhat intermediate between the B and C rings, 

562 strengthening the idea of a dependence on the optical depth. 



563 Finally, we conclude this section by remarking that the opposition effect is very diverse 

564 in Saturn's rings, and could be the consequence of different properties of the surface ring 

565 particles in various ring regions. Some general trends can be underlined, as we see in the 

566 next section. 

567 In a first step, we check if some correlations exist between the morphological parameters 

568 which depend on the same portion of the curve (e.g. A and HWHM for the surge) and 

569 also if parameters describing different parts of the phase curve can be correlated (e.g. A, 

570 HWHM and S for, respectively, the surge and the linear part). This is the purpose of the 

571 next section. 
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572 3.1.2 Cross comparisons between morphological parameters 

573 In order to constrain the morphology of the surge, we correlate the amplitude A with the 

574 angular width HWHM for the different main rings. 

575 We fit A and HWHM with a linear function and found correlation coefficients reported in 

576 Table 2. 



577 [Table 2] 



In Table 2 the Cassini Division, the amplitude A has the steepest function of the HWHM. 
A linear fit gives a slope of about 1.2 with a correlation coefficient of 47 %. This mean 
that narrow surges have a low amplitude and inversely wide surges have large amplitude. 
In the C ring, values of A and HWHM are generally greater than in the Cassini Division 
(but there are some exceptions). For the amplitude, this qualitati ve difference betwee n the 



583 Cassini Division an d the C ring is confirmed by previous study of 



results of 



Poulet et al. f2002h and 



French et al. f2007l). For the a ngular width, although our tren d of HWHM agrees w ith the 



Poulet et al. f2002h . this is 



French et al. (20071 ). it contradicts previous work of[ 
certainly due to the overestimation of HWHM when the smallest phase angles are missing 
(section 2.3.4 and Fig. 6b). 

For the C ring, we find a slope for A=f(HWHM) of about 0.9 with a good correlation 
coefficient of 79 %. 

The A ring shows a similar variation of A=f(HWHM) as for the C ring (1.0) with a 
correlation coefficient of 56 %. 

Finally, A and HWHM in the B ring have values smaller than in the faint rings (C ring 
and Cassini Division). Data points of A and HWHM for the B ring are concentrated 
in a similar range as the A ring one, however with a much shallower slope (0.6 with a 
correlation coefficient weakly reliable of 29 %, see Table 2). This means that the shape 
of opposition phase curves in the B ring may have various angular width with an almost 
constant value of the amplitude. 
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59B To conclude, the amplitude of the surge seems correlated with the angular width, at least 

599 for the C ring, A ring and Cassini Division whereas the amplitude is independent of 

600 HWHM for the B ring. The slope of A=f(HWHM) seems to decrease from the Cassini 

601 Division to the B ring, passing by the C and A rings, suggesting that the slope is a 

602 decreasing function of the optical depth. 

603 Whereas the slope S on the one hand, A and HWHM on the other hand are thought 

604 to be related to different portion of the phase curve (linear part and surge respectively), 

605 it is interesting to note that they are somewhat correlated. We simply note that S is a 

606 decreasing linear function the angular width and is also correlated with the amplitude. 

607 Slopes and correlation coefficients of S=f(A) and S=f(HWHM) are reported in Table 2. 

608 This coul d be due to the fact tha t we derive our slope from 1.5° to 25° whereas analytical 



609 model of 



Shkuratov et al. (19991 ). for example, describes the shado w hiding effect as a 



610 slope which fits the phase curve from 4.5° to larger phase angles, see (jPoulet et al.. 20021 ). 



611 3.1.3 Regional behavior 



612 We now look for regional behavior inside each ring (Fig. 8). The Fig. 8 displays A, HWHM 

613 and S as a function of the distance from Saturn by introd ucing a ring t ype nomenclature 



Cooke (19911 ) and decomposed 



614 based on the regional behavior of the C ring, well studied byL 

615 into three ring types : inner ring, background and plateaux. We have modified and extended 

616 this nomenclature to five classes of ring features, then applicable to the entire main rings 

617 system : 



6ib(1) inner regions characterized by low optical depth in all the rings (for example, the dark 



619 bands in the Cassini Division, see (IFlynn and Cuzzi. 19891 )) 



5(2) background are morphological smooth regions witho ut abrupt variation of r 



621 (3) bright regions [plateaux or plateaus in the C ring (IHolberg et al.. 19821 ). density and 



622 bending waves in the A ring located by 



Esposito et al. (1983| )) are the regions in each 
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623 ring with the highest optical depth, 



624(4) ringlets according to 
625 region or a gap, 



Holberg et aL (19821 ) are thinner ring embedded in a less dense 



626 ( 5) outer regions (for examp 



e the so-called ramp for the C ring and the Cassini Division 



627 (jCuzzi and Estrada. 19981 )) mark the transition at the boundaries of each ring. 



[Fig. 8] 



629 In Fig. 8a and b, the amplitude A and the HWHM vary smoothly across the main ring 

630 system (C to A rings), with only little scattering around the main trend (appart from 

631 the Cassini Division), illustrating the stability of the linear- by-part model for comparing 

632 multiple observations in different ring regions. From the C ring to the middle of the 

633 B ring, a smooth decrease is observed. No sharp transition is observed between the C and 

634 B rings. The outer regions of the C ring which are rich in gaps, plateaux and ringlets, 

635 have a somewhat larger value of amplitude. 

636 From the middle of the B ring to the outer of the A ring, both A and HWHM increase 

637 again. The Cassini Division presents (1) larger values of A and HWHM as in the C ring 
63B and (2) strongly dispersed values which may be indeed real because no image artifact is 

639 visible in this Division. 

640 The slope has a significantly different behavior (Fig. 8c) because strong jumps are observed 

641 at the boundaries of each ring. This reinforces differences of behavior of the surge and the 

642 linear part of the phase curve with the distance from Saturn. 



t, the behav ior of A, HWHM ans S did not show significant variations with 



643 As a resu 

644 ring type. ICooke (199ll ) noticed single scattering albedo of the C ring was dependent on 

645 the ring type classification. Maybe that the absence of correlation between morphological 

646 parameters and the ring type classification implies that the morphological parameters are 

647 independent of the single scattering albedo. S ince the single scattering albedo is correlated 



648 with the optical depth, (ISpilker et al.. 20051 ). we try now to correlate the morphological 
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649 parameters with the optical depth in the next section. 



650 3.1.4 Optical depth variations of the morphological parameters 

651 In order to quantify differences in terms of morphological shape in our 211 phase curves 

652 in CLEAR filters (see section 3.1.1), we use the three parameters A, HWHM and S. They 

653 are represented in the Fig. 9 as a function of the normal optical depth of the rings. 

654 [Fig. 9] 

655 Because the overall correlations of A, and S HWHM with r are not clear, although they 

656 seem to lead to an negative correlation for A and HWHM (Fig. 9a,b) and a positive 

657 correlation for S (Fig. 9c), we calculated with a linear fit the correlation coefficients for 
65B A=f(r), HWHM=f(r) and S=f(r) for each main ring, see Table 3. 

659 [Table 3] 

660 The amplitude A of the surge (Fig. 9a) is correlated with the optical depth of the rings. The 

661 following trends may be noted : first, the amplitudes of low optical depth (r<0.5, typically 

662 the C ring and the Cassini Division) are positively correlated with the optical depth (Ta- 

663 ble3) and second, the amplitude at high optical depth (t>1, typically the A and B rings) 

664 are negatively correlated with the optical depth (Table 3). The second trend was quite clear 

665 by eye, whereas the first one is more difficult to distinguish in Fig. 9a. This is due to the 

666 fact that the C ring and the Cassini Division have an optical depth restricted in the range 

667 0.01<r<0.5, which is quite compressed in the scale from to 2.5 of Fig. 9. The same be- 

668 havior (increasing of amplitude with increasing albed o, or r, the two values are correla ted 



Dovle et al.. 1989 : 



Belskaya and Shevchenko (20001 ) for 



Dones et al.. 19931 ) was noticed by[ 

670 the amplitude of asteroids' phase curves : for albedo<0.3 an increase of the amplitude with 

671 increasing albedo was noticed. 

672 We note also that for intermediate and high optical depth (0.5<r<2.5, typically the A 
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673 and B rings) the amplitude has a much smaller scattering and finally, the correlation co- 

674 efficient are larger (-74% for the B ring). 

675 For the A ring, it is interesting to note that regions of lower optical depth (0.3<r<0.5) 

676 connect well with data-points in the C ring both in terms in mean value and scattering 

677 (Fig. 9a). A good continuity with the B ring is also observed (0.7<r<l.l). 

67B The angular half width of the peak at half maximum (Fig. 9b) show similar correlation 

679 woth the optical depth than the amplitde of the surge : increase of HWHM when the op- 

680 tical depth increases for the C ring and the Cassini Division and decrease of HWHM 

681 when the optical depth increases for the A and B rings. The scattering of HWHM 

682 (HWHM~ 0.3 ± 0.2) for the low optical depth regions behaves similarly as for the ampli- 

683 tude, and lead to small correlation coefficients (19%, see Table 3). As for the amplitude, 

684 the scattering of HWHM is narrow for intermediate optical depth (0.5<r<1.6). Again, the 

685 behavior of the A ring is clearly intermediate between the C and B rings. To summarize, 

686 the behavior of HWHM is a decreasing function of increasing optical depth, with impor- 

687 tant scattering at low r that could be understood as an opposite behavior of HWHM 

688 (HWHM is an increasing function of increasing optical depth). 

689 The general trend for the slope of the linear regime (Fig. 9c) is a strong increase with in- 

690 creasing optical depth, with a uniform scattering and with central value well represented 

691 by S~ 0.07r^/^. This is the first time that a correlation is established between the slope S 

692 and the optical depth. Previous morphological study on asteroids' solar phase curves 

693 showed an exponential correlation between the slo pe of the phase function (the so- called 

694 phase coefficient (3 in mag.deg"^) and the albedo flBelskaya and Shevchenko. 2000l ). Our 

695 trend for the slope (the increase of S in wqP Aeg~^ with increasing r) is thus consistent 

696 with the slope of asteroids' phase functions (decrease of (3 in mag.deg"^ with increasing 

697 albedo, recall that the magnitude M is inversely proporti onal to I /F, the so-called geo- 



metric albedo : I /F = 10 o-^*^, see 



Domingue et al.. 1995 



so decrease of (3 in mag.deg ^ 



699 with increasing albedo leads to increase of /5 in I/F.deg ^ with increasing albedo), since we 
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ass ume that the optic a 



by 



Dovle et ah (19m : 



depth is positively correlated with the albedo, as already noticed 



Pones et al. (19931 ). However, as for A and HWHM, low optical 
depth regions (as the C ring and the Cassini Division) have a distinct trend than the trend 
of the A and B rings. 



704 In conclusion, whereas the slope S has a strong tendency with the optical depth (the 

705 average of the absolute correlation coefficients of S=f(r, see values of Table 3) is 53%), 

706 the first two parameters A and HWHM have a soft tendency with the optical depth (the 

707 average of the absolute correlation coefficients for A=f(r) is 44% and for HWHM=f(r) 
70B is 32%). Consequently, Fig. 9 which presents A, HWHM and S according to the optical 
709 depth, yields the following trends : 



710 (1) The morphological parameters of the surge (A and HWHM) have both similar behavior 

711 with r in the C ring and the Cassini Division. Firstly, a positive correlation of A with 

712 r and of HWHM with r, with a strong scattering and become almost independent 

713 of the optical depth for r>0.5. Secondly, we note a negative correlation of S with r. 

714 These trends are reversed (negative correlations of A and HWHM with r and a positive 

715 correlation of S with r) for the moderate and high optical depth regions (typically the 

716 A and B rings) ; 

717(2) the morphological parameter of the linear regime (S) is strongly positively correlated 

718 with the optical depth : negatively correlated in the C ring and the Cassini Division 

719 (where r<0.5) and positively correlated in the A and B rings (where r>0.5). 



720 We conclude that the trends of all the morphological parameters are linked to the optical 

721 depth. Then the drastic differences between, on the one hand, the amplitude and the 

722 angular width of the surge and, on the other hand, the slope of the linear part across 

723 the main ring system suggests that these characteristics originate from different physical 

724 mechanisms, as predicted by physical models (see section 2.3.5). 
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725 3.2 The opposition effect in COLOR filters 



726 3.2.1 Overview of the COLOR phase curves 

727 Fig. 10 details examples of phase curves obtained for the C, B, A rings and the Cassini 
72B Division. A striking result is that the COLOR phase curves' peaks are much narrower 

729 than the CLEAR phase curves' peaks. Indeed, this is due to the fact that WAC images 

730 (exclusively used for the CLEAR phase curves) do not have the more peaked part of the 

731 phase function (a>0.01°) whereas NAC images have it (a;>0.001°, see Fig. 4). To check 

732 this, we plot in the first column (labeled a) of Fig. 10 the phase function obtained with 

733 the camera WAC and in the second (labeled b) we plot the fuller phase function obtained 

734 with all the NAC and WAC images combined. Indeed, Fig. 10 demonstrates that the value 

735 of WqP{q) when a tend towards 0° is greater in the panel (b) than in the panel (a). 

736 [Fig. 10] 

737 The fact that WAC images do not include the peakest part of the opposition surge explains 
73B also the difference between the peak's intensity of curves in Fig. 3 in CLEAR filters and 

739 COLOR filters because when we processed images in CLEAR filters, there was no NAC 

740 images. This explains also the shift in I/F at the minimum phase angle on NAC images 

741 and WAC images of the same filter (Fig. 3). We observe that the shift in I/F between 

742 NAC and WAC image is minimum in the blue filter (Fig. 3), which implies that the full 

743 surge is contained in WAC images for this wavelength. 

744 In general, we noticed that the general shape of the curve is similar to that obtained 

745 previously : the C ring still has a fairly broad peak with a large amplitude (first pair of 

746 graphs at the top of Fig. 10). The B ring exhibits also narrow peaks as those of CLEAR 

747 filters comparatively to the other main rings (the second pair of graphs of Fig. 10). The 
74B Cassini division (third pair of graphs of Fig. 10) shows again a lot of scattering, which 
749 could be the consequence of the failure of the IChandrasekhar (1960l ) inversion when the 
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750 optical depth is not well known. Finally, the A ring, yet very dense, shows also scattering 



751 for WQP{a), which cou 

752 in Fig. 10. Indeed, the 



d be due to the gravi tational wakes, see the last pair of graphs 



Chandrasekhar (1960l ) inversion is very sen sitive to the optica. 



753 depth, and the wakes are known to modify locally the value of r f IColwell et al.. 2006 



Hedman et al. 2007bf ) 



755 3.2.2 Results in COLOR filters and comparisons with the morphological behaviors of the 

756 CLEAR filters 

757 We now consider the variation of morphological parameters A(A), HWHM(A) and S(A) in 
75B CLEAR and COLOR filters with the distance from Saturn (see Fig. 11). 

759 [Fig. 11] 

760 The amplitude of the opposition in the outer parts of the C and A rings has a much broader 

761 scattering in COLOR filters that in CLEAR filters (Fig. 11a). There is, unfortunately, not 

762 enough radial coverage to generalize this effect across the C ring. The Cassini Division 

763 also shows similar behavior for the dataset at high and low spectral resolution. In the 

764 B ring, where there is a good radial coverage in COLOR filters, amplitudes of CLEAR 

765 filters are much lower than the smallest amplitudes in COLOR (typically in the blue at 

766 A=0.451 /im), also noticed for the amplitudes in CLEAR and COLOR filters in the A 

767 ring. This could be due to the exclusive use of the images of the WAC to resolve the 
76B CLEAR opposition surge (section 3.2.1) 

769 The half-width at half maximum of the peak shows a good agreement between the data in 

770 CLEAR filters and COLOR filters (Fig. lib). The values are the same order of magnitude 

771 in all the rings and the regional effect is the same, except for the A ring where HWHM at 

772 high spectral resolution (COLOR filters) decreases when the distance at Saturn increases 

773 while at low spectral resolution (CLEAR filters), HWHM seems to start to grow. 

774 Finally, for the slope S of the linear part of the phase function, we represent a normalized 
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775 slope in order to compare data from the CLEA R filters with data from the COLOR filters. 



776 We follow the method of 



French et al. (20071 ) which divided the slope by the intensity of 



777 the background, thus the normalized slope is S/ Bi for the linear-by-part model and now 

77B has the deg~^ unit. In Fig. 11c we have for S similar behavior as a function from Saturn in 

779 the inner A ring and the Cassini Division for data in CLEAR filters and in COLOR filters. 

780 Finally, the B ring shows regional effects very different at high and low spectral resolution, 

781 especially in the middle of the B ring, at 107 000 km (the most optically thick region) 

782 where the CLEAR filters' slope is significantly overestimated compared to the slope in 

783 the blue filter (A=0.451 /im). In the same region, we also observed a strong scattering of 

784 HWHM (Fig. lib). 



785 3.2.3 Comparisons with multi-wavelength HST observations 



786 We now compa re the Cassini data (Fig. 11) with Earth-based data. For this purpose, the 



787 recent study of 



French et al. (20071 ) was chosen due to their small phase angles a>0.028° 



788 and one point below the n iinimum phase ang le corresponding to the angular size of the 



789 Sun. The phase curves of 



French et al. (20071 ) were obtained in I/F fo r the main rings 



Kaasalainen et al. (200ll ) which provided 



790 and adjusted with the linear-by-part model of 

791 morphological parameters A, HWHM and S for different wavelengths ranging from the 

792 ultraviolet to infrared (see their figure 7). 



French et al. f2007f ) 



793 We consider first the regional effects of morphological parameters derived by[ 

794 with the WFPC2/HST instrument. 

795 The dispersion observed by the HST for the three morphological parameters A, HWHM 

796 and S in the C ring and the Cassini Division is also very clear with ISS (Fig. llabc). 

797 However, it should be noted for the C ring scattering observed with our data (especially 

798 those in CLEAR filters) is highly localized (internal and external parts of the ring) and 

799 not present in the ce ntral regions, corresponding to the background according to the ring 



type classification of 



Cooke fl99l[ ). 
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In the B ring, [French et al. (20071 ) noticed strong variations of A and S in the inner edge 
which are absent for A and are present for S in the ISS data (Fig. 11a and c). For HWHM, 
HST did not have the raise of HWHM in the middle of the ring that we observed (Fig. lib). 
For the slope, we observe in both data set a decrease of S from the inner edge to the middle 
of the ring (Fig. 11c). 

le angular width H WHM are 



French et al. f2007f) obtain ed 



French et al. (20071 ) in 



806 In the A ring, the regional trends for the amplitude A and t 

807 similar for HST in ISS (Fig. 11a), but it is not the case for S. 

808 an almost constant value of S. The absolute value of the slope of 

809 all filters is roughly S~0.04 deg~^. With ISS/Cassini (Fig. 11c), we observe two distincts 

810 trends. First, S at high and low spectral resolution have similar values. Second, the regional 

811 effect observed is the same : there is a decreasing trend with increasing distance from Sat- 

812 urn. At the inner edge, S~0.025 deg~^ and at the outer edge S~0.015 deg~^. There is, 

813 however, a different order of magnitude obtained by ISS/Cassini and WFPC2/HST which 

814 could be due to lack of data at large phase angles for Earth-based observation. Indeed, 

815 the factor 2 between the slope values is well explained by the overestimation of the slope 

816 when the phase angle coverage stops at 6 degrees (section 2.3.4 and Fig. 6c). 



Now, we turn on the comparison on the variations of the morphological parameters with 
the w avelength, which lea ds to the following conclusions. 



First, 



of 



French et al. f2007l) found weak variations of A(A) in the B ring. Indeed, figure 3 



French et al. (20071 ) shows weak wavelength- variations of the amplitude in the B ring. 



and also in the A ring. Only one notable difference occurs in the ultraviolet. Overall, 
French et al. (200?! ) observed a decrease of A(A) from the ultraviolet to the green, and 



then an increase in the amplitude of the green to the infrared. This is not at all what 
is observed for the amplitudes of ISS where a notable increase from blue to green and a 
decrease in the infrared to the green (Fig. 11a). The values of A(A) for ISS/Cassini and 
WFPC2/HST a re anti-correlated. 



For HWHM(A), 



French et al. (20071 ) remarked a strong decrease of the half-width at half 
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82B maximum when the wavelength increases, and in all the rings. No break in decreasing 

829 HWHM(A) is observed by the HST. This decrease is very different of the behavior of 

830 HWHM(A) with ISS/Cassini (Fig. lib). Thus the variations of HWHM(A) presented here 

831 are also unprecedented. There is firstly not a decrease but an increase of HWHM with 

832 increasing wavelength. Second, the order of magnitude found is not the same. The angular 

833 half-width obtained with WFPC2/HST are generally between 0.05° and 0.2° while those 

834 of ISS are between 0.1° and 0.4°. This discrepancy could be the consequence of the relative 

835 thicknes s of our phase funct ions due to the use of several images at similar phase angles 



836 whereas 



French et al. (20071 ) did not have multiple recovered data because they obtained 



one point per phase angle. 



French et al. (20071 ) are consistent with those 



838 Finally, for S(A), the wavelength behavior of[ 

839 of ISS/Cassini (S decreases from the blue to the red, and slightly increases in the infrared), 

840 Fig. 11c. Strangely, this agreement shows a consistency of both data sets for variations 

841 of the slope with a wavelength whereas the regional effects of S observed by the two 

842 instruments are significantly different. 



843 3.2.4 Variations of the opposition effect with the wavelength and the optical depth 

844 We see that the morphological parameters depend on the optical depth (Fig. 9) and also 

845 on the wavelength (Fig. 12). 

846 [Fig. 12] 

847 It is appropriate now to quantify the variations of A(A), HWHM(A) and S(A) by fitting 
84B them with a linear model. We thus obtain for each COLOR phase curve of the rings two 

849 linear functions (one from 0.451 to 0.568 /im and another from 0.568 to 0.752 /im) which 

850 fit the behavior of A(A) and two linear functions HWHM(A), with the same wavelength 

851 boundaries. For the slope S(A), two linear functions are obtained (one from 0.451 yum to 

852 0.650 /im and from 0.650 /xm to 0.752 /im), see Fig. 13 for an example. 
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853 [Fig. 13] 

854 For this study, is only kept the slope of each linear function (that we called the steepness), 

855 that we correlated with the optical depth of the rings (the linear functions and their 

856 correlation coefficients are given in Tables 4 and 5). 

857 [Table 4, Table 5] 

858 The slopes of A(A) and HWHM(A) show both an increase with the optical depth from the 

859 blue to the green (Table 4). 

860 From the green to the infrared, we note a decrease of A(A) with r, but not for HWHM(A) 

861 which can increase or decrease with the optical depth. 

862 However, the strongest trends are found for S(A) which lead to an decrease with r of S(A) 

863 from the blue to the red and to an increase of S(A) from the red to the infrared. Values 

864 are given in Table 5. 

865 It is interesting to note the similar wavelength trends for A and HWHM, which are singu- 

866 larly distinct of the wavelength trends of S. This confirms the fact that the morphological 

867 parameters originate from different physical mechanisms, as predicted by the physical 
863 models (section 2.3.5). 

869 3.3 General trends of the opposition phase curves' morphology with the ISS data 

870 To conclude the section 3, the ISS data set provides several trends of the opposition effect 

871 in Saturn's rings which concern the : 

872(1) regional behavior of A, HWHM and S across the main ring system : indeed, 

873 the classification of ring type features defined in section 3.1.3 shows that A and HWHM 

874 vary continuously across the ring system (thus quite independently of the ring type 

875 classification), whereas S varies abruptly at the boundaries of each ring (Fig. 8) ; 

876 ( 2) morphology of the rings' phase curves : in section 3.1.4, strong dependences on 
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1 the optical depth of the rings of A, HWHM and S were estabhshed : anti-correlation of 
s A and HWHM with r and correlation of S with r (Fig. 9). A high scattering at low r 
) and high r is noticed for HWHM in CLEAR and COLOR phase curves ; 
)(3) correlation between the morphological parameters : we showed in the sec- 
L tion 3.1.2 that the parameters of the surge A and HWHM are correlated (Table 2), 
I whereas there is only a weak dependence of HWHM and S on the one hand, and A and 
? S on the other hand (Table 2) ; 

i(4) variations of the morphological parameters in CLEAR and COLOR filters : 

> we remarked in section 3.2.2 that the CLEAR amplitude is smaller than the COLOR 
5 ones, because the CLEAR data set did not includes NAC images which generally capture 
r th e highest part of the surge (Fig. 11). We showed also that the discrepancy in the slope 



i of [French et al. (20071 ) is due to their phase angle coverages which did not cover the 
) larger phase angles. As a consequence, their slopes are overestimated, as predicted by 
) our study on the influence of incomplete data (section 2.3.4 and Fig. 6c) ; 
[(5) variations of the morphological parameters with the wavelength and the 

> optical depth : A(A) increases from the blue to the green (0.451 /im-0.568 /im) and 
i decreases from the green to the infrared (0.568 yum-752 /im). Fig. 12. This increase is 
1 reinforced when the optical depth increases (Table 4) and is marked by a maximum value 

of A(A) in the green regardless of the rings. HWHM(A) also increases from the blue to 

> the green, but no general trend is noticed from the green to the infrared (Table 4). S(A) 
decreases from the blue to the red (0.451 /im-650 /im), then rises from the red to the 

i near-infrared (0.650 yum-752 /^m). Table 5. The wavelength- variation of S is then quite 
) distinct to that of the surge parameters (A and HWHM), then suggesting that they 
) may originate from different physical mechanisms. 
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901 4 Discussion 



902 ^.i A : a combination of the coherent backs cattering with the shadow hiding ? 



903 We showed in section 3.1.4 that the amplitude of the surge had a spe cific behavior ac- 



904 cording to the optical depth. With the model of 



Shkuratov et al. (19991 ) , an estimation of 



905 the amplitude is p ossible but refers only to the coherent backscattering enhancement. 



906 iPoulet et al. (20021 ) derived the following expression for the amplitude of the surge : 

907 y4 ~ 1 + e^'^/^/2 where d is the effective radius of grains and L is the free mean path of 

908 photons in the regolithic medium. It turns out that the amplitude could not be greater 

909 than 1.5, which contradicts our morphological results (Fig. 11a). This variation might be 

910 due to the shadow hiding effect which is not taken into account in this computation of 



the amplitude. 



912 4-2 Variations of S with the optical depth 



913 In shadow hiding numerical simulations of 



Stankevich et al. (19991 ) that take into account 



914 the optical depth and the filling factor, the absolute slope of the linear part has roughly 

915 the same value when the optical depth is greater than 1, whatever the filling factor value. 

916 This theoretical prediction is now confirmed by our results : we observed a saturation of 

917 S at r>l (Fig. 9c). 

91B Also, we observe with the ISS data (Fig. 11c) a wavelength -dependence of the slope S. 



919 This behavior was predicted by the model of 



Hapke (20021 ) because S depends on the 



920 extinction coefficient Qcxt(A, r) where r is the mean radius of particles (see section 2.3.4). 
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921 4-3 Linking photometric behaviors with dynamical concepts 



922 The present study, thanks to the quahty of its data (radial and angular resolution), allows 

923 us to highlight several observational facts never reported in the history of the observation 

924 of the opposition effect in the Saturn's system. The slope, the angular width and the 

925 amplitude of the rings' opposition phase curves are clearly correlated with the optical 

926 depth of the rings (Fig. 9). Whereas a physical description of such a dependence would 

927 need a full new physical model, we provide here some arguments explaining how the 
92B optical depth may indeed influence these three parameters. 

929 The optical depth is both a measure of the local volume filling factor of material, and 

930 the local collisional activity. Indeed, basic analytical computation shows that the filling 

931 factor of ring particle is proportional to t/H (with r standing for the optical depth and H 

932 standing for the vertical height of the rings). So one may expect regions of higher optical 

933 depth to have a much higher filling factor of particles. This is also in agreement with local 



934 simulation of ring dynamics flWisdom and Tremaine. 19881 ). that shows that the vertical 



935 width of material increases with decrease optical depth, because of the lower efficiency of 

936 collisional damping. So one may expect the filling factor of particles to be an increasing 

937 function of the optical depth. 



93B However, in the c ase the slope of the phase function for p 



939 numerical models (jlrvine. 1966 



Kawata and Irvine. 1974 



lase angles >1°, analyt ical and 



Stankevich et al.. 19991 ) predict 



940 that S should depend on the particle filling factor D and the vertical extension of the 

941 layer of particles in such way that the steeper is the slope, the lower is the filling fac- 

942 tor. Our morphological trends (Fig. 9c) seem unambiguously to contradict the theoretical 

943 trends of these models because our highest slopes are found in the high optical depth 



Salo and Karialainen. 2003 


), this implies that 


SD^l respectively for 


Stankevich et al., 1999 



assumptions of diluted lay ers (D<0.3 and 



Kawata and Irvine. 19741 ) are not suited 
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947 to the Saturn's rinffs. 

94B As a consequence, the shadow hiding observed in the Saturn's ri ngs may result of intra- 



949 particle shadow hiding, as stated by 



Buratti and Veverka (19851 ). Another possibihty is 



950 that the shadow hidin g that operates in the rings are mainta ined in a "regime of trans- 



951 parency" of particles (ILumme et al.. 1987 



Irvine et al.. 1988 



that dominates the shad- 



952 owing effect due to packing density. 



953 We now turn to the case of the HWHM and the amplitude. Whereas there is still a de 



3ate 



on what determine their value f Mishchenko. 1992 



Mishchenko and Dlugach. 1992a 



Shkuratov et al.. 199 



Hapke. 20021 ). authors agree to link them the coherent-backs cattering process, w hich may 



be controlled by the regolith at the surface of ring particles (jPoulet et al.. 20021 ). Like for 
planetary surfaces, the regolith is expected to be the result of space fractionation and 
erosional processes at the surface of the ring-particles, due (in particular) to the on-going 
collisional activity inside rings. Numerical studies of the dynamics of the ring-particles 
have shown that the optical depth is a key parameter controlling the collisional activity 



of rings. On the one hand, the number of col. 
to r (in the regime of low optical depth , see 



isions per orbit per particle i s proportional 



Wisdom and Tremaine. 19881 ). on the other 



hand, the random velocity in a ring of thickness H is about H xQ (with Q standing for the 
local keplerian frequency). Since if is a decreasing function of r, thus impact velocities 
are lower in regions of high optical depth. In short, particles in low-optical depth regions 
may suffer rare but violent collisions, conversely, in high-optical depth regions particles 
suffer frequent but gentle collisions. 

This may explain qualitatively why the HWHM and amplitude have different behavior in 
the data (Fig. 8 and 9). However, impact velocities have a lower bound ~ 2r x (with r 
standing for the particle's radius) due to the keplerian shear across the diameter of a par- 
ticle. This "shear dominated limit" is reached when the optical depth is high, typically for 
r>l. In such high filling factor regime the dynamics of collisions is entirely controlled by 
the keplerian shear rather than by the random impact velocities. This may qualitatively 
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974 explain why values of HWHM and Amplitude seem constant for r>l : in this regime, the 

975 collisional activity being about independent of optical depth, the physical properties of 

976 the regolith may be about constant which is indeed observed. 

977 In conclusion, in the absence of self-consistent physical model of the opposition effect, 
97B these qualitative arguments show that there are good reasons to believe that the optical 

979 depth is a key factor determining the opposition effect in the ring through two different 

980 mechanisms : 

981 (1) the optical depth may influence the absolute slope S, assuming that shadow hiding is 

982 the preponderant mechanism at phase angles a>l° 

983 ( 2) the optical depth may controls the HWHM and the amplitude (at phase angles a<l°) 

984 if the structure of the regolith is influenced by the collisional activity. 



985 4-4 Comparison of the opposition effect in optical light and infrared light 



986 A comparison of the solar phase curves of ISS/Cassini and the thermal phase curves of 

987 CIRS / Cassini yields to the following trends : 



988 (1) lAltobelli et al. (20071 ) found a prominent opposition surge for the thermal phase curves 

989 of the plateaux, well fitted by a logarithmic model. This is also the case of the solar 

990 phase curves of the p lateaux observed by ISS (Fig. 8a) ; 



991 (2) lAltobelli et al. (2007f ) found that the thermal phase curves of the background (regions 

992 in the close environment of the plateaux) do not have opposition surge, whereas back- 

993 ground has an opposition surge in the solar phase curves of ISS (Fig. 8a) ; 



994 A priori, the emitted phase curves may not reflect the coherent backscatter effect (lAltobelli et al.. 20071 ) 



995 because interferences of photons did no t act on heat, and thus on infrared light. However 



996 a pure shadow hiding model such as (ILumme and Bowell. 198ll ) fails to reproduce the 



997 CIRS opposition surge of the plateaux (in general the shadow hiding models did not 
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produce high surges, 



Stankevich et al.. 19991 ) . This could be the proof that the shadow 



999 hiding cannot produce solely the opposition surge in emitted light and that the coherent 

1000 backscatter could act on the shadow hiding mechanism, by r nultiplying th e single scat- 



Hapke (2002^. 



1001 tered light component at small phase angles, as underlined by 

1002 Interestingly, no similar surge were observed in the background (lAltobelh et al.. 20071 ) 

1003 whereas both plateaus and background have an opposition surge in the solar phase curves 

1004 of ISS. Because the background regions are more dim and reflective than the plateaus 



(background have smaller optical d epth and hig. 



they reflect more than they absorb. 



ler albedo than plateaus which means 



Cooke. 199X1 ). this could explain why these regions 



1007 did not have an opposition peak in emitted light. 



1008 5 Conclusion 



1009 We report here the main conclusions of this morphological study on the Saturn's rings 

1010 opposition effect seen by ISS/Cassini : 

1011 (1) The amplitude A and the half angular width HWHM of the opposition surge decrease 

1012 with increasing optical depth r. A and HWH M may reflect both coherent backscat- 



French et al. (20071 ) the morphological 



1013 tering and shadow hiding, because according to 

1014 parameters of the surge are greater than their coherent backscatter counterparts ; 
1015(2) All the morphological parameters are linked together. We find correlations between A 

1016 and HWHM, between HWHM and S, and also between A and S, which imply that 

1017 S could be more or less affected by the coherent backscattering. This could be due 

1018 to the fact that we d erive our slope from 1.5° to 25° whereas analytical model of 



Shkuratov et al. (19991 ) describes the shadow hi ding effect as a slo pe which fits the 



phase curve from 4.5° to larger phase angles, see (jPoulet et al.. 20021 ). 



1021(3) The absolute slope S of the linear part of the phase function increases with increasing 
1022 optical depth r (for optically thick rings) and shows distinct trends to the morpholog- 
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ical surge parameters, which implies that this parameter is not tota lly affected by the 



coherent backscattering. As (ILumme et al.. 1987 



Irvine et al.. 19881 ). we think that the 



i Saturn's rings could be in a regime of transparency of particles because the effect of 

5 packing density (decrease of S with increasing packing density) is not that expected for 

r the Saturn's rings (increase of S with optical depth, and if we assume that the optical 

1 depth and the packing density are correlated, increase of S with increasing packing 
J density) ; 

)(4) r-dependence with the morphological parameters strengthen our assumptions saying 

L that environmental effects are the key element determining the opposition effect because 

> the optical depth is a direct measure of the collisional and dynamical activity in the 
J surrounding of particles and is highly correlated with A, HWHM and S ; 

i(5) Comparisons of ISS/Cassini solar phase curves and CIRS/Cassini thermal phase curves 

> in the C ring show that the C ring's plateaus can have a strong opposition surge both in 
) solar and thermal phase curves whereas the C ring's background has a strong opposition 
r surge in the solar phase curve and no opposition surge in the thermal phase curve. 

J (6) Wavelength variations of the amplitude A of the surge show a maximum in all the rings 

) at A=0.568 fim. The increase of A from 0.451 to 0.568 fim and the decrease of A from 

) 0.568 to 0.752 fim are reinforced with increasing r ; 

l(7) Wavelength variations of HWHM of the surge show also a maximum at A=0.568 /im but 

2 it is not systematic, HWHM can also increase from 0.568 fim to 0.752 fim. Moreover, 
J there is no specific wavelength variations of HWHM with the optical depth ; 

i(8) Wavelength variations of the slope S of the linear part imply that the shadow hiding 

> depend on the wave length, may be via the particle's scattering cross-section of recent 



1046 model (IHapke. 20021 ). The decrease of S from 0.451 to 0.650 fim and the increase of S 



1047 from 0.650 to 0.752 fim are reinforced with increasing r. 



1048 The goal of this first paper was not to derive and quantify directly the physical properties 

1049 obtained from the models. First, because there is a large set of models and it seemed more 
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1050 convenient to separate the morphological models to the more physical and sophisticated 

1051 ones. Second, because recent physical models did not implement only the opposition ef- 

1052 feet but the main photometric effects which occur in the full phase function, from to 

1053 180 degrees. Consequently, more investigations will be provided for this purpose by using 

1054 full phase curves and photometric analytical models in the second paper. In the future, we 

1055 hope to have the linear degree of polarization at the opposition to obtain more constrains 

1056 on the ring particle' textures. 
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Table 1 

Main observational parameters of each sequence of images for each geometry of observation 
(i=arccos(/i) and e=arccos(/io))- CLEAR (WAC) filters designate the broadband filters in the 
optical domain (the central wavelength is A^'^'^ = 0.611/um). COLOR (NAC) filters designate 
blue, green, red and near infrared filters (at central wavelengths of X^^'^' = 0.440//m and A{^^° = 
0.451/im; A^.^^^ = 0.568/im; X^^^ = 0.650^m; Af*^ = 0.752^m) and COLOR (WAC) filters 
designate blue, green, red and near infrared filters (at central wavelengths of A^'^° = 0.460/^m; 
A^Ao = o.567^m; A*-^^ = 0.649/zm; A^^^ = 0.742^m). 



date 


Nb im 


i 


e 


Radial res. 


Azim. res. 


Filters (Camera) 






n 


(°) 


(km.pix "'") 


(km.pix "'" ) 




June 7 2005 


12 


111.5 


111.9 


44.0 


115.1 


CLEAR (WAC) 


June 26 2005 


66 


111.3 


111.3 


30.1 


70.0 


CLEAR (WAC) 


July 23 2006 


48 


73.1 


73.4 


13.4 


40.7 


CLEAR (WAC) 


May 20 2005 


57 


111.6 


111.6 


4.6 


11.5 


COLOR (NAC) 


May 20 2005 


59 


111.5 


111.9 


44.0 


115.1 


COLOR (WAC) 


Dec 31 2006 


12 


104.4 


108.7 


38.4 


104.8 


COLOR (WAC) 


Feh 20 2007 


20 


103.7 


122.7 


66.3 


72.5 


COLOR (WAC) 


Apr 25 2007 


16 


102.7 


110.7 


44.4 


92.6 


COLOR (WAC) 
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Table 2 

Results of linear fits (function and correlation coefficient) obtained for A=f(HWHM), S=f(A) 
and S=f(HWHM) for each Saturn's main ring and using CLEAR phase curves. 





A=f(HWHM) 


S=f(A) 


S=f(HWHM) 


function 


correl. 


function 


correl. 


function 


correl. 


Cass. Div. 


0.91 + 1.72r 


59 % 


0.04 -0.03r 


-68 % 


0.049 -O.OSr 


-35 % 


C ring 


1.20 + 0.91r 


79 % 


0.04 -0.06r 


-46 % 


0.025 -0.04t 


-47 % 


A ring 


1.12 + l.OSr 


56 % 


0.10 -0.14r 


-84 % 


0.079 -0.24t 


-81 % 


B ring 


1.18 + 0.58r 


31 % 


0.15 -0.33r 


-79 % 


0.141 -0.81T 


-65 % 
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Table 3 

Results of linear fits obtained (function and correlation coefficient) for A=f(T), HWHM=f(r) 



and S=f(r) for each Saturn's main ring and using CLEAR phase curves. 





A=f(T) 


HWHM=f(r) 


S=f(r) 


function 


correl. 


function 


correl. 


function 


correl. 


Cass. Div. 


1.22 + 0.937 


41 % 


0.20 + 0.175r 


19 % 


0.040 - 0.030T 


-42 % 


C ring 


1.50 + 0.024 


25 % 


0.33 + 0.077r 


19 % 


0.024 - 0.026T 


-68 % 


A ring 


1.50 - 0.175 


-36 % 


0.31 - 0.106r 


-49 % 


0.024 + 0.025T 


19 % 


B ring 


1.42 - 0.068 


-74 % 


0.24 - 0.016r 


-39 % 


0.035 + 0.031T 


82 % 
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Table 4 

Results of linear fits (function and correlation coefficient) obtained for the steepness of A(A)=f(r) 
and for the steepness of HWHM(A)=f(r) using COLOR phase curves of the A and B rings 
(r>0.5). 





Steepness of A(A)=f(r) 


Steepness of HWHM(A)=f(r) 




function 


correl. 


function 


correl. 


0.451 fiui <X< 0.568 i^im 


1.210 + 0.3r 


55 % 


0.14 + O.Ollr 


38 % 


0.568 fim <X< 0.752 j^im 


0.001 - 0.5r 


-51 % 


0.04 - 0.037T 


-46 % 
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Table 5 

Results of linear fits (function and correlation coefficient) obtained for the steepness of S(A)=f(r) 
and the steepness of S/Bi(A)=f(T) using COLOR phase curves of the A and B rings (r>0.5). 





Steepness of S(A)=f(r) 
(in wqP .deg~^ .iJ,m.~^) 


Steepness of S/Bi(A)=f(r) 
(in deg~^./im^^) 


function 


cor r el. 


function 


correl. 


0.451 <A< 0.650 
0.650 fim <A< 0.752 fim 


-0.21 - 0.043r 
0.04 + O.OlSr 


-73 % 
67 % 


-0.04 - 0.039T 
0.02 + 0.054T 


-49 % 
44 % 
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Fig. 1. The opposition effect in the B ring. A typical image of the 26 June sequence captured by 
the Wide Angle Camera (W1498453136.IMG). Concentric circles correspond to identical phase 
curves computed in the image (isophase). Numbers in the concentric circles represent the solar 
phase angles (in degrees). 



59 



Fig. 2. The opposition effect in the C ring. A typical image of the 20 may sequence captured 
by the Narrow Angle Camera (N1595278165.IMG). The contrast is enhanced to make more 
visible the opposition spot in the C ring. Concentric circles correspond to identical phase curves 
computed in the image (isophase). Numbers in the concentric circles represent the solar phase 
angles (in degrees). 
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Fig. 3. Radial location of the opposition spot in the images taken in CLEAR and COLOR filters 
(filled symbols represent WAC images and empty symbols represent NAC images). We give in 
the y-axis the normalized brightness I/F of the minimum phase angle in the image, the x-axis 
is the corresponding distance from Saturn of this point. The vertical dotted lines correspond to 
ring boundaries. The optical depth of PPS/Voyager is plotted as a radial reference. 
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Fig. 4. Extracted COLOR phase curve from NAC images of the May 20 sequence in A, B, 
C rings and Cassini Division in linear scale (a) and logarithmic scale (b). The vertical dotted 



line in (a) and (b) corresponds to the angular size of the Sun. The solid curves in (a) and (b) 

3). 



correspond to the logarithmic model of 



Bobrov ri97C 



62 



(a) Typical B ring phase curve (R = 110 340 km) (b) Typical C ring phase curve (R = 96 496 knn) 




5 10 15 20 5 10 15 20 

Phase angle a (°) Phase angle a 



Fig. 5. Typical phase curves of the B ring (a) an d C ring (b) in CLEAR filters chosen for 



Kaasalainen et al. (200ll ) when the phase 



testing the stability of linear-exponential model of 
angle coverage is incomplete. Circles and squares correspond to the cutoff of the incomplete 
ph ase function and the soli d curves correspond to the initial fit of the linear-exponential model 
of iKaasalainen et al. (200 ll ). The vertical dotted lines correspond to boundaries of the surge 



domain and the linear part domain that we have tested (see Fig. 6). 
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Fig. 6. Deviation of the hnear-exponential model of 
logical parameters A (a), HWHM (b) and S (c) when the phase angle coverage is incomplete. 
Circles and squares correspond to the cutoff of the incomplete phase functions of respectively 
the B ring and C ring typical phase curves (Fig. 5). 
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Fig. 7. Representative CLEAR phase curves for the main rings with a zooi n on the surge (a) 



and the full phase curve (b), fitt ed respectively with the lo garithmic model of lBobrov (197d ) (a) 



and the linear- by-part model of iLumme and Irvine (19761 ) (b) 
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Fig. 8. Regional behavior of the morphological parameters : (a) the amplitude A, (b) the angular 
width H WHM in degrees and (c) the absolute slope S in rooP-deg"^ from the Linear-by-part 



model of iLumme and Irvine (19761 ) for CLEAR phase curves using the ring type classification. 
The vertical dotted lines correspond to ring boundaries. The optical depth of PPS/Voyager is 
plotted in (c) as a radial reference. 
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Fig. 9. Morphological pa rameters of CLEAR phase curves from the Linear-by-part model of 



Lumme and Irvine (19761 ) : Amplitude A (a), Angular width HWHM (b) and absolute slope 



S (c) in roo-P.deg ^. 
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Fig. 10. Representative COLOR phase curves for the main rings with a zoom on the surge with 
WAC images (a) and the full phase curve with NA G and WAG images (b). Fuh phase curves of 



(b) are fitted with the linear-exponential model of iKaasalainen et al. (20011 ) 
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Fig. 11. Regional behavior of morphological parameters : (a) the amplitude A, (b) the angular 
width HWHM and (c) un itless absolute slope S/Bi in deg~^ from the Linear-by-part model of 



Lumme and Irvine (19761 ) using CLEAR and COLOR phase curves. The vertical dotted lines 
correspond to ring boundaries. The optical depth of PPS/Voyager is plotted in (c) as a radial 
reference. 
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Fig. 12. Morphological parameters of CLEAR and COLOR phase curves from the Linear-by-part 



model of 



Lumme and Irvine (19761 ) : (a) the amplitude A, (b) the angular width HWHM and 



(c) the unitless absolute slope S/Bi in deg ^. 
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Fig. 13. Variations of the morphological parameters A(A), HWHM(A), the absolute slope S(A) 
and the unitless absolute slope S/Bi{X) for two typical regions of the A ring (a) and in the B 
ring (b). Dotted lines correspond to linear fits obtained in the spectral ranges. The slopes of 
these linear functions are called "steepness" and are correlated with the rings' optical depth in 
tables 4 and 5. 
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Table 1: Outputs of the logarithmic model of lBobrov (1970l ) for the CLEAR phase curves 
representing each ring types [inner corresponds to inner regions characterized by low 
optical depth, background are morphological smooth regions without abrupt variation of 
optical depth, bright corresponds to bright regions that have the highest optical depth 
in each ring, ringlet corresponds to a thinner ring embedded in a less dense region or 
a gap, and outer corresponds to outer regions that mark the transition at each ring 
boundary) of each main ring. Horizontal lines correspond to the ring boundaries. We give 
for each ring type, the normal optical depth of the Voyager/PPS instrument and the 
corresponding radius (the distance from Saturn's center). frooP^i/F is a conversion factor 
which corresponds to the mean level of I/F curves over the mean level of wqP curves. 



Rad. (km) 


''"pps 


Ring type 


ao 




frooP^I/F 


73271.7 


0.005 


ringlet 


2.023 


-0.6457 


0.00001 


74996.2 





058 


bright 


1 


105 


-0 


2097 


0.0223 


75356.9 





027 


inner 


1 


331 


-0 


2857 


0.1255 


75665.0 





032 


inner 


1 


072 


-0 


2535 


0.1232 


76199.8 





132 


bright 


1 


128 


-0 


1899 


0.1255 


76671.2 





031 


inner 


1 


188 


-0 


2561 


0.1253 


77075.4 





119 


bright 


1 


263 


-0 


2106 


0.0643 


77862.3 





731 


ringlet 





568 


-0 


1074 


0.0488 


78273.0 





079 


background 


1 


242 


-0 


2130 


0.0501 


78889.1 





074 


background 


1 


182 


-0 


2030 


0.0545 


79238.2 





331 


bright 





807 


-0 


1237 


0.0310 


79956.9 





096 


background 


1 


217 


-0 


1985 


0.1189 


80757.8 





109 


background 


1 


213 


-0 


1962 


0.1259 


table continues on next page... 
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Rad. (km) 


''"pps 


Ring type 






frooP-^I/F 


81661.4 


0.117 


background 


1.186 


-0.1845 


0.0466 


82031.0 


0.202 


bright 


0.969 


-0.1540 


0.0539 


82914.0 


0.133 


background 


1.242 


-0.1907 


0.0621 


84125.6 


0.102 


background 


1.352 


-0.2052 


0.0557 


84844.4 


0.425 


bright 


1.707 


-0.2325 


0.0418 


85234.6 


0.099 


background 


1.290 


-0.1953 


0.0432 


85706.9 


0.256 


bright 


1.335 


-0.1978 


0.1254 


85953.3 


0.227 


bright 


1.098 


-0.1625 


0.1252 


86158.7 


0.075 


background 


1.506 


-0.2179 


0.1249 


86503.2 


0.396 


bright 


1.853 


-0.2430 


0.1250 


86877.4 


0.066 


background 


1.404 


-0.1984 


0.0522 


87189.3 


0.153 


bright 


0.892 


-0.1213 


0.0430 


87312.8 


0.163 


bright 


1.019 


-0.1440 


0.1140 


87506.6 


1.011 


ringlet 


0.645 


-0.0566 


0.1163 


88451.6 


0.239 


bright 


1.869 


-0.2257 


0.1241 


88725.6 


0.156 


ringlet 


1.087 


-0.1390 


0.1223 


89233.7 


0.248 


bright 


1.695 


-0.2038 


0.1221 


89547.1 


0.045 


background 


1.995 


-0.2505 


0.1232 


89851.1 


0.307 


bright 


1.970 


-0.2265 


0.1233 


90019.4 


0.073 


background 


2.361 


-0.2384 


0.0687 


90163.1 


0.713 


ringlet 


1.373 


-0.1719 


0.0633 


90509.7 


0.355 


bright 


2.310 


-0.2514 


0.0895 


90685.9 


0.076 


outer 


2.454 


-0.2532 


0.1144 


90929.2 


0.099 


outer 


1.981 


-0.2161 


0.0985 


91237.3 


0.137 


outer 


2.060 


-0.2253 


0.1128 
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Rad. (km) 


''"pps 


Ring type 






frooP-^I/F 


91788.6 


0.170 


outer 


2.276 


-0.2447 


0.0857 


92390.9 


1.629 


bright 


3.548 


-0.3668 


0.0952 


93781.3 


1.359 


inner 


3.398 


-0.3418 


0.1254 


94201.6 


0.935 


inner 


3.377 


-0.3479 


0.1256 


94751.3 


0.706 


inner 


3.010 


-0.3528 


0.1253 


95107.0 


0.732 


ringlet 


2.984 


-0.3582 


0.1249 


95527.3 


0.936 


inner 


3.421 


-0.3460 


0.1252 


96659.0 


1.340 


inner 


3.637 


-0.3472 


0.0542 


96885.4 


0.817 


ringlet 


3.213 


-0.3492 


0.0608 


97661.4 


0.976 


inner 


3.413 


-0.3394 


0.1265 


98146.4 


0.872 


inner 


3.497 


-0.3387 


0.1264 


99213.5 


1.458 


background 


3.721 


-0.3418 


0.1259 


100475. 


2.057 


bright 


3.882 


-0.3486 


0.1281 


101671. 


1.188 


ringlet 


3.698 


-0.3279 


0.0057 


101800. 


2.104 


bright 


3.939 


-0.3401 


0.1259 


103223. 


1.266 


background 


3.975 


-0.3634 


0.1236 


105422. 


2.070 


bright 


4.149 


-0.3446 


0.0818 


106327. 


2.069 


background 


4.139 


-0.3577 


0.0477 


107847. 


2.116 


bright 


4.302 


-0.3634 


0.0679 


109367. 


2.053 


bright 


4.428 


-0.4181 


0.1249 


110789. 


1.815 


ringlet 


4.250 


-0.4146 


0.1251 


111953. 


1.407 


background 


4.292 


-0.4617 


0.0766 


112309. 


1.910 


outer 


4.412 


-0.4370 


0.1249 


113053. 


1.768 


outer 


4.441 


-0.4350 


0.1253 


113441. 


1.744 


outer 


4.393 


-0.4334 


0.1263 


table continues on next page... 



74 



Rad. (km) 


''"pps 


Ring type 






frooP-»I/F 


113667. 


1.358 


background 


4.306 


-0.4675 


0.1269 


113796. 


1.563 


outer 


4.394 


-0.4632 


0.1261 


114411. 


1.746 


outer 


4.344 


-0.4653 


0.1262 


114637. 


1.440 


background 


4.342 


-0.4818 


0.1266 


115122. 


2.047 


bright 


4.506 


-0.5140 


0.1268 


116254. 


2.105 


bright 


4.547 


-0.4917 


0.1249 


117803. 





698 


ringlet 





722 


-0 


2112 


0.1265 


117910. 





087 


ringlet 


1 


506 


-0 


0942 


0.1266 


117983. 





113 


inner 


2 


119 


-0 


2572 


0.1272 


118084. 





119 


inner 


1 


883 


-0 


2250 


0.1269 


118168. 





150 


bright 


1 


329 


-0 


2249 


0.1288 


118241. 





099 


ringlet 


1 


324 


-0 


1754 


0.0204 


118365. 





080 


inner 


2 


172 


-0 


1818 


0.0222 


118482. 





091 


inner 


1 


960 


-0 


2528 


0.0464 


118668. 





079 


inner 


1 


897 


-0 


2426 


0.0225 


118836. 





083 


inner 


1 


718 


-0 


2306 


0.0475 


119061. 





082 


bright 


1 


878 


-0 


2338 


0.0435 


119145. 





089 


inner 


1 


958 


-0 


2494 


0.1003 


119229. 





104 


bright 


1 


640 


-0 


2175 


0.0508 


119285. 





097 


bright 


1 


824 


-0 


2418 


0.0561 


119476. 





029 


inner 


2 


604 


-0 


3045 


0.0626 


119644. 





038 


inner 


1 


986 


-0 


2493 


0.0527 


119768. 





032 


inner 


2 


186 


-0 


2975 


0.0639 


120060. 





356 


ringlet 


1 


509 


-0 


2411 


0.0518 


120116. 





161 


background 


1 


939 


-0 


2629 


0.0560 
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Rad. (km) 


''"pps 


Ring type 






frooP-^I/F 


120279. 


0.134 


ringlet 


1.600 


-0.2347 


0.1158 


120335. 


0.077 


background 


1.980 


-0.2603 


0.0448 


120408. 


0.064 


background 


2.168 


-0.2863 


0.0423 


120565. 


0.367 


background 


2.328 


-0.3048 


0.0370 


120638. 


0.385 


background 


2.574 


-0.3377 


0.0051 


120711. 


0.384 


background 


2.625 


-0.3319 


0.0326 


120773. 


0.502 


background 


2.621 


-0.3319 


0.0465 


120918. 


0.089 


outer 


2.637 


-0.3581 


0.0952 


121031. 


0.116 


outer 


2.318 


-0.2747 


0.1058 


121272. 


0.140 


outer 


2.427 


-0.2736 


0.1093 


121626. 


0.156 


outer 


2.685 


-0.2951 


0.0563 


121901. 


0.180 


outer 


2.786 


-0.3296 


0.0790 


122097. 


0.458 


inner 


2.887 


-0.3562 


0.1254 


122269. 


0.623 


inner 


3.366 


-0.4469 


0.1253 


122553. 


1.092 


bright 


3.862 


-0.4958 


0.1247 


123040. 


0.864 


inner 


3.696 


-0.4431 


0.1254 


123249. 


0.681 


inner 


3.478 


-0.4386 


0.1235 


123676. 


1.231 


bright 


3.911 


-0.4554 


0.1254 


123848. 


0.951 


bright 


3.731 


-0.4466 


0.1251 


124252. 


0.878 


bright 


3.576 


-0.4149 


0.1238 


124409. 


0.690 


bright 


3.451 


-0.4045 


0.1246 


124659. 


0.574 


background 


3.305 


-0.3920 


0.1245 


125367. 


0.699 


bright 


3.075 


-0.3468 


0.1256 


125951. 


0.590 


bright 


3.212 


-0.3655 


0.1256 


126619. 


0.498 


background 


3.074 


-0.3872 


0.1256 
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Rad. (km) 


''"pps 


Ring type 






frooP-»I/F 


127655. 


0.495 


background 


2.929 


-0.3742 


0.1257 


128977. 


0.473 


background 


2.694 


-0.3465 


0.1256 


130037. 


0.476 


background 


2.568 


-0.3344 


0.1254 


130786. 


0.661 


bright 


2.637 


-0.3575 


0.1254 


131122. 


0.465 


background 


2.556 


-0.3406 


0.1254 


131818. 


0.803 


bright 


2.764 


-0.3932 


0.1253 


132372. 


0.735 


bright 


2.478 


-0.3015 


0.1253 


133105. 


0.463 


background 


2.438 


-0.3082 


0.1253 


133513. 


0.012 


ringlet 


3.354 


-1.0378 


0.1253 


133574. 


0.027 


ringlet 


1.028 


-0.3205 


0.1251 


133820. 


0.509 


outer 


2.196 


-0.1977 


0.1250 


133917. 


0.534 


outer 


2.375 


-0.2581 


0.1250 


134280. 


0.959 


bright 


2.390 


-0.2465 


0.1249 


134493. 


0.575 


outer 


2.435 


-0.2647 


0.1249 


134737. 


0.585 


outer 


2.440 


-0.2657 


0.1247 


134965. 


0.558 


outer 


2.429 


-0.2550 


0.1246 


135166. 


0.548 


outer 


2.340 


-0.2079 


0.1245 


135595. 


0.614 


outer 


2.282 


-0.1600 


0.1247 


135808. 


0.596 


outer 


2.232 


-0.1423 


0.1246 


135913. 


0.523 


outer 


2.263 


-0.1494 


0.0234 


136070. 


0.586 


outer 


2.194 


-0.1160 


0.0526 


136288. 


0.602 


outer 


2.154 


-0.0992 


0.1254 


136579. 


0.627 


outer 


1.795 


-0.0080 


0.1216 


136665. 


0.590 


outer 


1.963 


-0.0159 


0.1206 


136736. 


0.824 


outer 


1.714 


-0.0040 


0.1198 
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Rad. (km) 


''"pps 


Ring type 








140338. 


0.065 


ringlet 


0.240 


-0.0090 


0.1214 
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